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CHAPTER  1 


INTRODUCTION 


A  knowledge  of  snow  growth  is  important  not  only  in  the 
understanding  of  wintertime  precipitation  itself.  but  also  in 
other  kinds  of  precipitation  when  the  Wegener-Bergeron 
process  is  important.  Besides  the  understanding  of 
precipitation  processes.  the  investigation  of  snow  growth 
mechanisms  is  also  essential  to  artificial  cloud  seeding  and 
to  small  scale  parameterization  in  large-scale  general 
circulation  models  and  in  numerical  weather  prediction 
models. 

The  study  of  snow  growth  can  be  by  different  means,  such 
as  laboratory  investigat ion*  numerical  modelling.  in  situ 
observations  and  mathematical  analysis.  Each  of  these 
approaches  has  its  advantages  and  limitations.  Laboratory 
investigations  on  snow  growth  have  the  advantage  of  being 
able  to  study  it  in  a  controlled  environment.  The 
temperature,  humidity  and  other  factors  can  be  under  strict 
control  of  the  experimenter.  However,  the  conditions  inside 
a  laboratory  are  rather  ideal  and  can  be  quite  different  from 


15 


the  actual  atmospheric  environment.  In  addition*  because  of 
the  difficulty  of  performing  laboratory  experiments*  the 
physical  processes  of  snow  growth  can  be  studied  only  one  at 
a  time.  Laboratory  investigations  on  secondary  ice  crystal 
production  by  crystal-crystal  collision  has  been  performed  by 
Vardiman  (1978)  who  stresses  the  possiblity  of  collisional 
breakup  in  explaining  the  presence  of  the  high  density  of  ice 
crystals  in  stratiform  clouds  with  embedded  convection.  The 
drawback  of  this  experiment  is  that  ice  crystals  are  casued 
to  impact  on  hard  surface  which  is  quite  different  from  the 
effect  of  cry stal-crystal  collision.  Blanchard  and  Spencer 
(1970)  perform  experiments  on  the  the  binary  interactions  of 
raindrops.  They  establish  the  importance  of  drop  breakup  in 
determining  an  exponential  drop  size  distribution.  Mossop 
and  Mallet  (1974)#  Mallet  and  Mossop  (1974)  and  Mossop  (1976) 
report  laboratory  experiments  on  secondary  ice  crystal 
production  due  to  rime-splintering.  They  conclude  that  the 
rate  of  splinter  production  is  a  function  of  the 
concentration  of  supercooled  liquid  water  droplets  larger 
than  24  microns  between  -8*C  and  -3*C. 

Numerical  modelling  has  the  advantage  of  being  able  to 
incorporate  many  physical  processes  of  snow  growth  plus  the 
interaction  between  microphysics  and  dynamics  in  a  single 
study.  The  drawback  of  numerical  modelling  is  that  physical 
insights  can  easily  be  lost  in  a  large  model.  Moreover*  the 


computation  time  required  for  a  full  model  can  be  enormous. 
Young  (1975)  performs  numerical  model  of  raindrop  evolution 
with  the  physical  processes  of  condensation,  coalescence  and 
breakup.  It  is  found  that  the  effect  of  collisional  breakup 
dominates  over  the  effect  of  spontaneous  breakup  and  produces 
an  exponential  raindrop  size  distribution.  The  results  are 
insensitive  to  the  collision  efficiency  and  the  steady-state 
spectrum  is  insensitive  to  the  number  of  breakup  fragment. 
Gillespie  and  List  (1976)  study  the  evolution  of  raindrop 
size  distributions  in  steady  state  rainshafts  using  a 
numerical  stochastic  coalescence-breakup  model.  The  results 
conclude  that  there  is  only  one  equilibrium  distribution 
which  depends  on  the  rainfall  rate  for  each  Marshal 1-Palmer 
distribution.  Leighton  (1980)  uses  a  numerical  model  to 
study  the  depostion  and  aggregation  processes  of  snowflakes 
and  concludes  that  the  assumption  of  exponential  snow  size 
distribution  is  valid  and  he  also  confirms  the  analytical 
results  of  Passarelli  (1978b). 


In  situ  observations  of  snow  growth  became  feasible 
since  the  development  of  laser  imaging  probes.  This 
technique  of  studying  snow  growth  has  one  advantage  which  no 
other  techniques  have  in  that  what  is  observed  is  also  what 
actually  happens  in  the  atmosphere.  However.  the  data 
obtained.  such  as  the  snow  size  distribution  (e.  g.  Lo  and 
Passarelli.  1982a)  and  the  two  dimensional  shapes  of  snow 
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crystals  <e. g.  Dyer  and  Glass.  1982).  are  the  results  of 
snow  growth,  while  the  physical  processes  of  snow  growth  have 
to  be  inferred  from  these  data.  The  changes  of  snow-size 
spectra  with  height  have  been  performed  by  Passarelli  (1978a. 
c>  and  by  Houze  et  al.  (1979).  The  findings  show  that  snow 
size  spectra  change  rather  systematical ly  with  height.  A  way 
to  interpret  the  change  of  snow-size  spectra  will  be 
presented  in  chapter  2.  Gordon  et  al  (1982)  uses  a  technique 
similar  to  Lo  and  Passarelli  (1982a).  They  conclude  that 
snow  growth  goes  through  four  stages  where  different  physical 
processes  predominate,  namely  the  vapor  deposition  stage,  the 
aggregation  stage,  the  secondary  ice  crystal  production  stage 
and  the  an  aggregation  stage. 

Analytical  studies  of  snow  growth  provides  the  best 
insight  concerning  the  physical  mechanism  of  snow  growth. 
But.  in  order  to  simplify  the  mathematics  so  that  Analytical 
solutions  can  be  obtained,  many  assumptions  have  to  be  made 
in  the  set  up  of  the  problem.  Also,  in  order  to  keep  the 
mathematical  problem  tractable.  only  a  few  of  the  physical 
processes  can  be  treated  at  one  time.  Srivastava  (1971) 
studies  the  effect  of  coalescence  and  breakup  in  raindrop 
growth  and  concludes  that  the  binary  interactions  give  rise 
to  an  equilibrium  size  distribution  slope.  Passarelli 
(1978b)  calculated  the  evolution  of  snow-size  distribution  in 
a  winter-time  stratiform  snow  situation  using  an  analytical 


model  with  vapor  deposition  and  aggregation.  He  concludes 
that  the  equilibrium  snow  size  spectra  owe  their  existence  to 
the  counteracting  effects  of  deposition  and  aggregation 
growth. 

Previously*  the  main  thrust  of  the  study  of  the 
microphysical  processes  of  precipitation  growth  has  been  on 
raindrops.  Few  investigations  have  been  made  to  study  the 
snow  growth  processes  more  comprehensively.  The  study 
presented  here  is  an  attempt  to  give  a  comprehensive  picture 
of  the  snow  growth  processes  once  the  snow  crystals  have  been 
initiated.  The  methodology  used  is  basically  a  combination 
of  all  four  techniques*  with  particular  emphases  on  in  situ 
aircraft  observations  and  mathematical  analyses. 

Chapter  2  presents  the  in  situ  observational  data  and 
their  analyses.  The  observations  are  made  by  using  a  single 
aircraft.  This  new  flight  plan  is  devised  so  as  to  permit 
one  to  observe  the  height  evolution  of  snow-size  spectra  in  a 
reference  frame  where  the  effects  of  horizontal  gradients  and 
temporal  changes  in  the  atmosphere  are  minimized  (Lo  and 
Passarelli*  1982a*  see  appendix  1).  This  flight  plan*  termed 
the  advecting  spiral  descent  (ASD>«  requires  an  aircraft  to 
start  aloft  in  a  mesoscale  precipitation  area  (e. g.  a  snow 
shaft*  Marshall*  1993)  and  spiral  downward  in  a  constant  bank 
angle*  descending  at  approx imately  the  mean  fallspeed  of 
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snow.  The  loops  of  the  spiral  drift  with  the  wind.  The 
analysis  is  performed  by  averaging  spectra  over  a  complete 
loop  of  the  spiral*  serving  to  average  any  horizontal 
inhomogeneities • 

The  snow-size  spectra  are  first  analysed  using  the 
conventional  negative  exponential  size  distribution.  From 
the  behavior  of  the  size  spectra*  the  physical  processes  of 
snow  growth  can  be  inferred.  From  these  analyses>  the 
horizontal  gradients  of  microphysical  structure  are  also 
studied.  The  data  are  then  analysed  using  a  negative 
exponential  power  law  snow-size  distribution.  This  new 
technique  of  describing  the  size  distribution  provides 
certain  new  insights.  Finally*  these  data  are  compared  with 
other  flight  data. 

Chapter  3  presents  theoretical  models  of  snow  growth. 
The  physical  processes  represented  in  the  models  include 
vapor  deposition*  aggregation.  collisional  breakup  and 
rime-splintering.  The  basic  equation  makes  use  of  the 
stochastic  collection  formulation.  Snow  size  spectra  are 
represented  either  in  a  negative  exponential  distribution  or 
in  a  negative  exponential*  power  law  distribution  (Lo  and 
Passarelli*  1982b >.  From  the  basic  equation*  moment 
conservation  equations  are  derived.  The  equations  are 
simplified  analytically  as  far  as  possible  and  the  resulting 
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Chapter  2 

ASD  DATA  AND  ANALYSES 

2.  1  ASD  Data 

Advccting  spiral  descants  were  flown  with  an  Air  Force 
Geophysics  Laboratory  C-130  aircraft  that  is  equipped  with 
thermodynamic  and  cloud  physics  sensors  (Barnes  at. al. > 
1982).  The  thermodynamic  sensors  are  a  dew  point  hygrometer 
and  total  air  temperature  probe.  The  cloud  physics  sensors 
include  a  Ewer  probe  for  total  water  and  ice*  icing  rate 
detector*  Johnsons-Wi 1 1 iams  (J-W)  cloud  water  probe»  total 
water  content  indicator*  formvar  hydrometeor  replicator  and 
several  Particle  Measuring  Systems  (PMS)  probes>  including  an 
axially  scattering  spectrometer  probe  (ASSP>*  a  1-D  cloud 
probe<  a  1-D  precipitation  probe*  a  2-D  cloud  probe  and  a  2-D 
precipitation  probe  (Knol lenberg*  1970). 

In  this  study*  basically  only  the  data  from  the  1-D 
precipitation  probe  have  been  used  in  the  quantitative 
analyses.  This  is  a  laser  imaging  probe  that  counts  and 
sizes  particles  into  15  size  categories  ranging  from  300  to 
4500  microns.  The  data  from  the  1-D  cloud  probe  (which 
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measures  particles  in  the  size  range  of  20  to  300  microns) 
are  not  used  because  there  is  a  serious  discrepancy  in  the 
number  density  at  the  size  range  where  the  cloud  probe  and 
the  precipitaton  probe  overlap.  It  is  true  that  certain 
investigators  have  attempted  to  smooth  the  differences  but 
the  smoothing  techniques  tend  to  be  rather  artificial  (Belksy 
et. al. •  1981).  The  2-0  data  are  used  qualitatively  to 
determine  the  crystal  type.  No  quantitative  analyses  of  the 
2-0  data  have  been  done.  The  data  from  the  J-W  water  content 
probe<  icing  rate  detector  probe  and  ASSR  are  used  to 
determine  cloud  liquid  water  content.  However>  there  is  a 
problem  inherent  with  the  ASSP  in  snow-storm  situations.  The 
ASSP  operates  on  a  light  scattering  principle  and  therefore 
can  easily  give  false  counts  from  the  light  scattered  or 
reflected  off  the  snow  particles. 

There  is  a  difficulty  in  using  the  airborne  temperature 
measurement  in  that  the  true  air  temperature  computed  is 
sometimes  lower  than  the  dew  point  temperature  by 
approximatedly  1* C.  In  plotting  the  soundings>  saturation  is 
assumed  whenever  this  happens. 

Extensive  analyses  from  three  flights  are  presented  with 
supp lementary  data  from  other  flights.  Two  of  the  three 
flights  were  made  off  the  coast  of  Washington  (25  and  26 
February  1980)  and  one  off  the  coast  of  New  Hampshire  (8 


March  1980).  One  ASO  was  performed  on  each  day  and  these 
will  be  referred  to  as  spirals  1>  2  and  3  respectively.  All 
three  spirals  took  place  in  winter  cyclonic  storms.  Spiral  1 
was  flown  just  ahead  of  an  occluded  front  (Fig.  2.  1)  and 
spiral  2  was  flown  approx imately  100  km  east  of  a 
low-pressure  center  (Fig. 2. 2).  Spiral  3  was  flown  in  the 
warm  sector  of  a  cyclonic  storm  (Fig.  2.  3).  The  soundings 
obtained  from  the  aircraft  data  (Fig.  2.  4)  indicate  that  the 
atmosphere  was  slightly  more  stable  than  moist  adiabatic  and 
was  saturated  or  near  saturated  in  all  cases.  The  height  and 
temperature  ranges  for  the  three  spirals  are  given  in  Table 
2.  1. 

Snow  size  spectra  averaged  over  a  complete  loop  are 
shown  in  Figs.  2.  S.  2.6  and  2.7  for  spirals  1.2  and  3 
respectively.  The  spectra  are  labeled  consecutively  from  the 
top  to  the  bottom  loops  of  the  spiral.  Because  the  spectra 
are  almost  negative  exponential,  the  plots  on  a  linear-log 
scale  result  in  straight  lines.  Therefore,  the  concept  of 
intercept  and  slope  can  be  applied  here.  The  spectra  from 
the  first  spiral  show  an  initial  monotonic  increase  in  the 
intercept  while  the  slope  remains  relatively  constant  through 
loop  12.  after  which  there  is  a  rapid  decrease  in  both  the 
slope  and  intercept  during  loops  13  and  14.  The  remainder  of 
the  loops  show  essentially  no  change.  The  vertical 
separation  between  successive  loops  is  200  m. 
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The  spectra  from  the  second  spiral  show  a  similar 
behavior.  The  intercept  shows  an  initial  montonic  increase 
while  the  slope  remains  approx imately  constant  from  loops  1 
through  7.  Then  both  the  intercept  and  the  slope  decrease 
from  loops  7  through  14.  The  spectra  of  loops  IS  to  19  are 
almost  constant.  The  vertical  separation  between  successive 
loops  is  also  200  m  for  this  case. 


The  spectra  from  the  third  spiral  indicate  a  somewhat 
different  behavior.  There  was  not  a  stage  where  the 
intercept  increases  while  the  slope  remains  constant.  Both 
the  intercept  and  the  slope  decrease  from  loop  1  to  loop  10. 
Then  the  slope  and  the  intercept  remain  more  or  less  constant 
for  loops  11  through  14.  The  vertical  separation  between 
successive  loops  is  about  100  m  for  this  case. 


The  data  from  the  J-W  cloud  water  probe  are  extremely 
noisy.  So<  the  occurrence  of  supercooled  liquid  water  is 
determined  using  the  Rosemont  icing  rate  detector  data  and 
the  ASSP  data.  Fig.  2.  8  shows  the  icing  rate  measurements  for 
the  spiral  on  25  February  1980.  The  icing  probe  is  a  rod 
protruding  outside  the  aircraft.  In  the  'dry'  state,  the  rod 
has  a  certain  natural  vibration  frequency.  When  there  is 
supercooled  liquid  water  in  the  atmosphere,  icing  would 
develop  on  the  rod,  thus  changing  the  mass  of  the  rod. 
Therefore,  the  vibration  frequency  of  the  rod  changes.  The 


sec.)  elapsed  after  the  beginning  of  flight. 


Fig. 2. 9: 


ASSP  particle  concentration  (with  smallest  size  category  ignored) 
for  25  February,  1980 
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liquid  water  content  is  proportional  to  the  change  in 
frequency  which  is  in  turn  proportional  to  the  derivative  of 
the  trace  shown  in  the  figures.  In  order  to  smooth  the  data, 
a  nine  minute  running  mean  technique  is  applied.  The  data 
are  plotted  against  the  time  elapsed  since  the  beginning  of 
the  spiral.  The  data  indicate  that  supercooled  liquid  water 
is  present  at  the  end  of  the  spiral  (^3500  sec.  after  the 
start  of  the  spiral).  The  time  when  the  snow  growth 
mechanism  changes  is  also  marked  on  the  graph.  The 
concentration  of  particles  detected  by  the  ASSP  (with  the 
smallest  size  category  ignored  and  also  smoothed  by  a  nine 
minute  running  mean)  (Fig.  2.  9)  indicates  huge  increases 
around  3300  sec.  This  again  suggests  the  existence  of 
supercooled  liquid  water  and/or  secondary  production  of  ice 
particles. 

The  data  of  the  icing  rate  detector  for  the  spiral  on  26 
February  1980  are  platted  in  Fig.  2.  10.  They  indicate  that 
there  is  liquid  water  at  1300  sec.  .  then  at  3300  sec..  then 
at  3900  sec.  and  at  4200  sec.  Finally,  at  4500  sec.  the 
icing  rate  detector  goes  through  a  complete  cycle.  These 
show  that  patches  of  supercooled  liquid  water  can  be  found 
throughout  most  of  the  flight  and  the  liquid  water  content  is 
even  higher  at  the  end.  The  ASSP  concentration  is  plotted  in 
Fig.  2.  11.  The  data  show  that  there  is  no  change  in 
concentration  until  3300  sec.  and  3900  sec.  However,  there 
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Fig. 2. 10:  Icing  rate  measurement  for  26  February,  1980  (Layout  similar 
to  Fig. 2. 8) 
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Fig. 2. 11:  *.SSP  particle  concentration  (with  smallest  size  category  ignored) 
for  26  February,  1980 
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i*  no  big  increase  in  concentration  at  4500  sec.  when  the 
icing  rate  detector  goes  through  a  complete  cycle.  The  icing 
rate  detector  is  located  about  3  m  from  the  ASSP.  which 
could  explain  the  apparent  incoherence.  The  overall 
concentration  for  this  spiral  is  lower  than  that  of  25 
February  1980. 

The  icing  rate  detector  data  and  ASSP  detected 
concentration  for  the  spiral  on  8  March  1980  are  plotted  in 
Figs.  2.  12  and  2.13.  From  the  graphs>  no  liquid  water  is 
evident  for  this  spiral.  In  summary,  supercooled  liquid 
water  is  present  but  rare  in  some  of  the  flights  but  its 
occurrence  is  not  correlated  with  any  particular  stage  of 
snow  growth. 

2.  2  Two-parameter  data  analyses 

Because  of  the  quasi-exponential  behavior  it  is  possible 
to  approximate  the  spectra  by  an  exponential  form  such  that 
(Qunn  and  Marshall.  1958) 

N(D)dD  -  N0e'XI>  dD  (2.  1) 

where  N(0)dD  is  the  particle  concentration  in  the  diameter 
interval  CD.  D+dDl.  N,  is  the  intercept  and  the 
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Fig. 2. 12:  Icing  rate  measurement  for  8  March,  1980  (layout  similar 
to  Fig. 2. 8) 
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Fig. 2. 13:  ASSP  particle  concentration  (with  smallest  size  category 
ignored)  for  8  March,  1980 


distribution  slope.  Here  D  is  the  actual  1-dimensional 
particle  size  measured  by  the  laser  imaging  probe.  Using  a 
least-squares  regression  technique!  Ne  and  X  for  each 
spectrum  can  be  found.  In  order  to  ensure  adequate 
statistics!  it  is  stipulated  that  for  any  size  category  to  be 
included  in  the  regression  calculation  the  total  number  of 
particles  counted  in  that  size  category  must  be  greater  than 
10  for  the  complete  loop.  Fig.  2.  14  shouts  the  data  from  the 
three  spirals  in  log,,  N.-log^/V  space  (hereafter  N0  -A  space). 
The  loop  numbers  are  indicated  next  to  the  data  paints  which 
are  connected  sequentially  by  a  line.  Environmental 
temperatures  are  indicated  for  places  where  significant 
changes  occur. 


This  type  of  spectral  display  is  convenient  because  a 
straight  line  in  space  corresponds  to  a  constant  moment 
of  an  exponential  distribution!  since  the  jth  moment  of  a 
spectrum  is 


M;  -  f  D>  N,  e'APdD  -  — - 


(2.  2) 


Hence  in  logarithmic  N,  -Aspacei  the  moment  is  constant 

along  any  straight  line  having  slope  (j+l>. 


The  behavior  of  the  spectra  in  Figs.  2.  5i  2.6  and  2.7  can 
be  discussed  in  terms  of  the  trajectory  of  the  spectral 
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•volution  in  N0~\  space  in  Fig.  2.  14.  For  spiral  1,  there  is 
a  gradual  increase  in  Ns  and  a  slight  decrease  in  /\  from 
loop  1  thorugh  loop  12  such  that  the  slope  of  the 
No-A  trajectory  is  negative!  implying  that  all  spectral 
moments  are  increasing.  At  loop  12  the  spectral  evolution 
changes  dramatically  and  both  N,  and  X  decrease  until  loop 
13.  After  that  N(  and  A  assume  approx  imately  constant 
values.  The  second  spiral  shows  a  similar  pattern  of  three 
stages  of  evolution.  Stage  1  (loops  1  through  7)  is 
character i zed  by  an  increase  in  N.  accompanied  by  relatively 
little  change  in  A  .  Stage  2  (loops  7  through  14>  is 
character i zed  by  a  rapid  decrease  in  both  N,  and  A  .  Stage  3 
is  marked  by  an  apparent  cessation  of  spectral  evolution. 
The  third  spiral  does  not  show  the  presence  of  stage  1.  Both 
N.  and  A  decrease  from  loops  1  through  10.  Then  loops  10 
through  14  indicate  that  the  spectral  evolution  has  come  to  a 
halt. 


The  N# - A  trajectories  for  the  second  stage  of  growth 
(both  intercept  and  slope  decreasing)  are  roughly  parallel 
for  all  three  spirals«  having  slopes  ranging  from  1.80  to 
1.93.  This  suggests  that  during  this  phase  of  spectral 
evolution!  the  sum  of  the  diameters  of  snow  particles  is  a 
conservative  property  of  the  d istr ibution.  These  three  cases 
also  suggest  that  the  distribution  slope!  has  a  minimum 
value  of  approximately  10  cm'1  .  which  character i zes  the  third 
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stag*  of  evolution. 


2.  3  Interpretation  of  two-parameter  data  analyses 

The  effect  of  depositan  growth  on  the  size  distribution 
is  character i zed  by  how  the  rate  of  change  of  the  particle 
diameter  depends  on  the  particle  diameter/  i.  e.  , 

JjL  oo  dS  (2.  3> 

dt 

i 

For  Of  larger  particles  grow  more  quickly  in  diameter 

f 

space  and  one  would  expect  the  distribution  slope  to 
decrease.  For  0#  depositional  growth  will  cause  the 

distribution  to  steepen.  For  snow#  the  particle  mass  varies 
approximately  with  the  square  of  the  particle  diameter  and 
the  rate  of  change  of  mass  is  directly  proportional  to  the 
particle  diameter  so  that  $2  0.  Hence/  to  a  first 
approximation/  all  sizes  will  grow  at  the  same  rate  in 
diameter  space  and  the  distribution  slope  will  remain 
constant.  The  intercept  will  increase  since  the  smaller/ 
more  numerous  particles  will  grow  into  larger  sizes.  In  the 
absence  of  a  source  of  small  particles/  a  lower  limit  to  the 
size  distribution  will  develop. 

The  first  stage  of  particle  growth  (observed  for  spirals 


1  and  2)  is  consistent  with  deposition  growth.  This  is  not 
to  say  that  other  physical  processes  are  not  occurring  but 
that  deposition  is  dominating.  Also  in  this  stage*  the 
increase  in  N#  and  slight  decrease  in  ^  indicate  an  increase 
in  the  total  concentration  of  particles.  This  is  probably 
due  to  the  growth  of  small  particles  into  sizes  detectable  by 
the  200-Y  probe.  The  source  of  these  small  particles  could 
be  either  nucleation  and/or  secondary  production. 

The  second  stage  of  growth*  where  both  and 
/S  decrease*  is  observed  for  all  three  spirals.  This  is 
character istic  of  aggregation  which  depletes  small  particles 
and  creates  large  ones.  The  ASD  data  for  spirals  1  and  2 
show  a  very  sudden  transition  from  the  dominance  of 
deposition  to  aggregation.  During  the  deposition  growth 
phase  the  particles  are  small  and  collisions  are  rare* 
eventually  deposition  produces  a  sufficient  number  of  large 
particles  for  aggregation  to  commence.  Once  started.  the 
large  particles  produced  by  aggregation  accelerate  the 
aggregation  process  and  rapidly  deplete  the  smaller 
particles.  This  accounts  for  the  sudden  transition  and  the 
very  rapid  evolution  of  the  size  distribution  after  the 
transition.  For  spirals  1  and  2  the  transition  occurs  at 
-13*C  which  is  in  the  dendritic  growth  temperature  range. 
This  implies  that  the  tendency  for  dendrites  to  form 
aggregates  <e. g. *  Jiusta  and  Weickmann.  1973)  may  play  a  role 


in  the  transition  to  aggregation  growth.  However  the  model 
calculations  described  by  Passarelli  (1978a,  b,c>  and  the 
computations  presented  in  the  next  section  do  not  require  a 
change  in  the  particle  geometry  or  collision  efficiency  to 
simulate  the  first  two  stages  of  growth. 

The  most  puzzling  behavior  is  the  apparent  sudden  end  to 
aggregation.  The  spectra  cease  evolution  when  the  slope 
reaches  approx imately  10  cm-1  .  This  suggests  that  the 
depletion  of  small  particles  is  balanced  by  a  production 
mechanism.  This  mechanism  can  be  either  primary  production 
or  secondary  production  or  both.  The  temperature  range  for 
this  stage  is  0aC  to  -10'C,  where  ice  nucleation  is  not 
predominant  (Flectcher,  1962)  and  so  the  probability  of 
primary  production  is  low.  Moreover  large  particles  produced 
by  aggregation  are  somehow  depleted.  This  suggests  that 
breakup  is  occurring.  In  veiw  of  the  presence  of  liquid 
water,  there  is  also  the  possibility  of  rime-splintering. 

The  fact  that  all  three  spirals  evolve  to  the  same  slope 
suggests  collisional  breakup.  This  hypothesis  is  drawn  from 
previous  work  on  drop  coalescence  and  breakup  which 
demonstrates  that  coalescence  and  collisional  breakup  leads 
to  equilibrium  distributions  which  have  the  same  slope, 
regardless  of  the  precipitation  rate  (Gillespie  and  List, 
1976  and  Srivastava,  1978).  The  latter  author  concludes  that 
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raindrop  size  distributions  tend  to  have  approx imate lg  a 
constant  slope  and  an  intercept  proportional  to  rainfall 
content.  This  is  in  agreement  with  Blanchard  and  Spencer's 
(1970)  observations  on  raindrop  size  distributions  in  which 
binarg  interaction  processes  would  produce  exponential 
distributions  with  a  constant  slope. 


2.4  Horizontal  gradients  in  the  microphgsical  structure 

In  order  to  examine  the  horizontal  gradient  in  the  size 
distribution*  the  mean  diameter  of  all  particles  measured  bg 
the  200— Y  spectra  as  a  function  of  aircraft  magnetic  heading 
is  plotted  in  Fig.  2.  15  for  the  first  and  second  spirals. 
Each  loop  is  represented  bg  a  separate  graph  and  the  loops 
are  stacked  verticallg  in  accordance  with  their  height.  The 
approximate  diameter  of  each  loop  is  6  km.  The  loop  in  which 
the  spectra  transformed  from  first  to  second  stage  of  growth 
(the  peak  in  the  N0-)t  trajectorg)  is  indicated  bg  a  star  in 
each  case.  Before  the  transition*  the  mean  diameter  is 
essentiallg  uniform  within  each  loop  and  graduallg  increases 
with  depth.  (Not  all  upper-level  loops  are  shown).  However* 
after  the  transition*  horizontal  inhomogeneities  develop  verg 
rapidlg.  The  features  are  correlated  from  one  loop  to  the 


next. 
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On*  interesting  feature  in  the  second  spiral  is  that  at 
loop  9/  the  average  particle  diameter  has  its  maximum  around 
320  magnetic  heading.  Then*  this  maximum  shifts  gradually 
towards  the  west  and  then  the  south.  This  perhaps  indicates 
that  the  wind  shear  is  turning  cyclonically  downward  (or 
anticyc lonical ly  upward)  and  the  snow  particles  are  being 
carried  around  with  it. 

To  examine  the  point-to-point  behavior  of  the  spectra* 
w*  arbitrarily  divide  each  loop  into  four  quadrants  bounded 
by  the  cardinal  directions.  Ideally  the  spectral  evolution 
can  be  studied  by  examining  the  height  evolution  of  spectra 
averaged  over  a  particular  quadrant  rather  than  over  an 
entire  loop.  Fig.  2.  16  shows  an  example  of  the  N,-X  evolution 
for  the  second  spiral  for  the  four  quadrants.  The  N, scale  is 
different  for  each  quadrant  in  order  to  separate  the  four 
quadrants.  Note  that  the  four  N, -A  trajectories  all  show  the 
same  general  features.  However*  the  west-north  quadrant 
starts  rapid  stage  2  growth  earlier  than  the  the  other 
quadrants  (*.  g.  >  examine  loops  9  and  10  in  Fig.  2.  IS). 

Fig.  2.  16  illustrates  that  even  if  we  examine  portions  of 
a  loop  the  spectral  evolution  is  coherent.  This  is 
consistent  with  the  previous  discussions.  Also*  the  rapid 
development  of  horizontal  gradients  of  mean  diameter  are 
apparently  related  to  the  fact  that  spectra  in  different 


regions  undergo  the  transition  from  deposition  to  aggregation 
at  different  heights. 


2.  5  Three-parameter  analyses 

Snow-size  spectra  have  conventionally  been  represented 
by  a  two-parameter  negative  exponential  distribution  (2.1). 
This  expression  has  the  advantage  of  being  mathematically 
simple.  However*  this  expression  is  not  a  good  description 
of  the  small  end  of  the  spectrum.  A  three-parameter  negative 
exponential-power  law  distribution 

T 

N(D) dD  *  N^(h)e  D  dD.  (2.4) 

where  N(D)dD  is  the  number  density  of  particles  with 
diameters  in  the  range  0  and  D+dD.  provides  another  degree  of 
freedom.  Hence*  it  can  provide  a  better  fit  to  the  observed 
snow-size  spectra*  even  at  the  small  end.  Such 
three-parameter  distr ibutions  have  been  used  by  Takeuchi 
(1978)  and  Uplinger  (1981)  to  describe  raindrop  size  spectra. 
Any  theoretical  analysis  using  a  three-parameter  size 
distribution  will  also  have  one  more  dimension  to  manipulate. 

Discretizing  (2.4)  into 
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I°g(0  Na  *  logBN4  -X4Di  log„e  +  <rlog(#  D i  +  error.  (2.5) 

Let  Ki  -  Da 

*  logw  Nx 
a  *  log„  N, 
b  -  X4log,0  e 
c  =  o~ 

«  *  number  of  size  categories 


Using  least-squares  fit  regressions*  a.  b  and  c  can  be  found 
from  the  spectra  to  be 


l  %  u  i  *  • 

~  ~  “  r 


(2.  6) 


b  =  F,  +  cF^ 


<2.  7) 


where 
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The  results  for  N^f ram  spiral  1  is  shown  in  Fig.  2.  17. 
N the  'intercept'  goes  through  three  distinct  stages  as 
snow  grows.  First#  N*.  increases  slowly.  Then  N^.  takes  a 
rapid  rise  and  a  rapid  fall.  Finally  remains  relatively 
constant.  These  three  stages  correspond  to  the  deposition# 
aggregation  and  collisional  breakup  stages  in  the 
two-parameter  analyses.  Fig.  2.  18  shows  the  change  of  N*  with 
height  for  the  second  spiral.  The  data  appear  to  be 
extremely  noisy  especially  at  the  beginning  of  the  spiral. 
This  is  probably  due  to  the  lower  particle  number  density  in 
that  portion  of  the  spiral.  Also  note  that  the  scale  in  this 
figure  is  different  from  that  for  spiral  1.  Still,  the 
general  behavior  of  N*  is  similar  to  that  of  the  first 
spiral.  first  increases#  then  it  decreases  and  finally 
takes  on  a  somewhat  constant  value. 

Fig.  2.  19  shows  the  evolution  of  with  height  for  the 
third  spiral.  The  result  for  this  spiral  is  rather  different 
from  those  of  the  previous  two  spirals.  As  recalled  from  the 
two-parameter  analyses  (Fig.  2.  14>#  this  is  the  spiral  in 


growth  are  marked  by  dashed  lines. 


-  26-80 


Fig. 2. 18:  N*  as  a  function  of  height,  26  February,  1980,  similar  to  Fig. 2. 17 
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which  there  is  no  stage  1  (where  NQ  increases  and  X  stags 
more  or  less  constant).  Except  for  the  third  and  fifth 
loops*  is  generally  decreasing.  Comparing  this  with  the 
two-parameter  model*  loops  1  thorugh  10  correspond  to  stage  2 
and  loops  10  through  14  correspond  to  stage  3.  That  means  in 
stage  2*  where  aggregation  is  the  dominant  effect,  small 
particles  are  being  depleted  and  N#  is  decreasing.  In  stage 
3.  where  aggregation  and  collisional  breakup  come  to  an 
equilibrium.  N^,  is  approximately  constant. 

Figs.  2. 20.  2.21  and  2.22  show  the  change  of  X* .  the 
exponent.  with  height  for  spirals  1.  2  and  3  respectively. 
X^decreaes  monotonical ly  as  the  snow  crystals  fall.  This  is 
similar  to  the  behavior  of  X  in  the  two-parameter  analyses. 
In  the  last  few  loops  of  each  of  the  spirals.  \  would  become 
more  or  less  constant.  However,  unlike  the  two-parameter 
analyses,  the  terminal  values  of  X*  are  different  for  each  of 
the  spirals.  This  is  probably  due  to  the  fact  that  X* is 
sensitive  to  the  value  of  T  .  As  will  be  seen  later.  <T  is 
quite  different  for  each  of  the  spirals. 

Fig.  2.  23  gives  the  change  of  <r  with  height  for  spiral  1. 
The  observed  values  indicate  that  in  the  early  part  of  the 
flight.  <r  is  negative  and  has  an  approximately  constant 
value.  Suddenly,  0"  changes  to  become  positive  and  then  it 
takes  on  a  more  or  less  constant  value.  The  place  where 
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cr  makes  the  sudden  change  from  negative  to  positive 
corresponds  to  the  middle  of  the  aggregation  stage.  A 
negative  <T  means  that  there  are  many  more  small  particles 
than  large  particles  and  a  positive  <r  means  that  the  number 
of  small  particles  is  limited  and  there  is  a  size  at  which 
the  particle  density  peaks  (Fig.  2.  24).  The  sudden  change 
from  negative  to  positive  implies  that  at  the  aggregation 
stage.  small  particles  are  depleted  very  rapidly. 

Fig.  2.  25  shows  the  change  of  <r  with  height  for  sprial  2. 
Just  as  the  case  of  N*  for  the  same  sprial.  the  values  of 
(Tare  very  noisy.  If  we  ignore  the  first  four  loops.  is 
initially  negative.  Then  it  increases  to  become  positive  and 
then  stays  apporximately  constant.  This  is  similar  to  that 
of  spiral  1. 

Fig.  2.  26  shows  <r  versus  height  for  spiral  3.  Here  the 
evolution  of  <f  is  vastly  different  from  the  previous  two 
cases.  Instead  of  having  an  initially  negative  value.  r  is 
initially  positive  and  then  it  decreases  into  the  negative 
from  loop  1  through  loop  10  (corresponding  to  stage  2  in  the 
two-parameter  analysis).  <r  remains  negative  from  loops  10 
through  14  (corresponding  to  the  third  stage  in  the 
two-parameter  analyses).  Physically*  this  means  that  the 
number  density  of  small  particles  increases  faster  than  that 
of  the  larger  particles.  This  could  imply  either  the 
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probability  of  collisional  breakup  is  high  or  the  number  of 
fragments  generated  during  collisional  breakup  is  extremely 
high  or  both.  This  question  will  be  investigated  in  the 
section  on  theoretical  analyses. 

2.6  Comparisons  with  other  studies 

Some  recent  observational  studies  of  snow-size  spectra 
have  shown  a  spectral  behavior  similar  to  the  ASD  results/ 
although  the  scatter  is  typically  much  greater.  This  is 
possibly  due  to  the  fact  that  the  data  are  not  taken  in  a 
manner  as  coherent  as  the  ASD.  Passarelli  ( 1978a# c)  employed 
similar  instrumentation  but  a  different  vertical  sampling 
scheme.  The  aircraft  was  placed  in  either  a  constant  ascent 
or  descent  while  flying  a  constant  heading.  Particle  size 
distributions  were  averaged  over  15  s  intervals.  The  spectra 
from  flights  on  6  and  10  March#  1975  are  shown  in  N.-X  space 
in  Fig.  2.  27.  While  spectra  obtained  via  the  ASD  technique 
show  a  very  systematic  behavior  with  height#  the  spectra 
obtained  on  these  two  days  do  not.  The  spectra  apparently 
lie  on  a  line  corresponding  to  the  second  stage*  but  the 
position  is  random.  The  spectra  for  each  day  are  very  well 
differentiated  in  space  which  probably  reflects  the  very 
different  environmental  conditions  on  the  two  days. 


o  6  Mor  1975  (Poseorelli  1978a  ) 

•  10  Mar  1975  (  "  "  ) 

*— *  26  Nov  1975  (  '•  1978c) 

A  22  jon  1975  (Houze  etal.1980) 

□  8  Dec  1976  (  "  «t  oi  1980 ) 


Related  observational  studies.  See  text  for  details 
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Passarelli  (1978c)  also  averaged  particle  size  spectra 
over  10  km  horizontal  passes  at  various  altitudes.  The  data 
for  a  flight  on  26  November  1975  are  also  shown  in  Fig.  2.  27. 
The  passes  are  spaced  at  600  m  and  timed  at  10  min  apart* 
which  is  approximately  the  time  required  for  snout  to  fall 
from  one  level  to  the  next.  The  results  shout  a  first  and 
second  stage.  Vertical  incidence  radar  measurements  on  this 
day  indicated  fairly  steady  precipitation  while  the  aircraft 
was  sampling. 

Houze  et  al.  *  (1979*  1980)  employed  a  similar  technique 
by  flying  level  passes  although  no  attempt  was  made  to  follow 
snow  from  level  to  level.  The  spectra  for  flights  on  22 
January  1976  and  8  December  1976  are  also  shown  in  Fig.  2.  27. 
The  data  are  more  scattered*  perhaps  because  transient 
environmental  conditions  are  manifested  over  the  long 
sampling  paths.  However*  the  observed  spectra  are  still 
within  the  values  from  the  other  cases.  Although  the  height 
for  each  spectrum  is  not  indicated*  Houze  et  al.  (1979)  show 
that  N0and  X  decrease  with  increasing  temperature.  Based  on 
the  general  trend  that  temperature  increases  with  decreasing 
altitude*  it  can  be  inferred  that  N. and  X  decrease  with 
decreasing  altitude*  which  is  in  agreement  with  the  second 
stage  concept. 


Oordon  et  al.  (1982)  flew  an  elongated  figure-eight 
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pattern  while  the  plane  was  descending  at  1  or  2  m/s. 
approximately  the  fallspeed  of  snow.  The  figure-eight  was 
oriented  normal  to  the  orientation  of  a  rain  band.  Their 
analyses  technique  is  similar  to  that  used  for  the  ASD.  The 
data  are  derived  from  a  PMS  2D-P  precipitation  probe*  but  it 
is  not  clear  how  the  2-dimensional  data  are  reduced  to  1 
dimension.  Their  results  for  15  February  1982  are  plotted  in 
Fig.  2. 28.  From  the  results  in  Na-)t  space*  their  findings  are 
similar  to  those  from  the  spirals.  However*  they  conclude 
that  there  are  four  stages  of  snow  growth.  The  first  stage, 
at  temperature  below  -22° C*  snow  growth  is  dominated  by 
nucleation  and  deposition*  which  basically  agrees  with  the 
findings  from  the  ASD  data.  The  second  stage,  between  -11°  C 
and  -22*C  is  where  aggregation  is  dominant*  which  again 
agrees  with  the  ASD  data.  The  third  stage*  from  -4°C  to 
-10°C«  is  where  secondary  ice  crystal  production.  and 
aggregation  are  dominant.  In  contrast  to  this  study,  these 
authors  conclude  that  secondary  ice  crystal  production  is  due 
to  a  rime-splintering  mechanism.  They  base  their  reasoning 
for  rime-splintering  on  the  concentration  of  particles 
(assumed  to  be  water  drops)  larger  than  24  microns  as  counted 
by  the  Forward  Scattering  Spectrometer  Probe  (FSSP).  As  has 
been  presented  earlier*  the  particle  counts  from  the  FSSP 
during  snow  conditions  are  rather  unreliable.  Gordon  et  al. 
did  not  use  other  instrments*  such  as  the  J-W  probe  or  icing 
rate  detector  to  double  check  the  presence  of  liquid  water. 
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So#  th«  occurrence  of  r irae-sp 1  inter ing  is  not  well 
substantiated.  The  fourth  stage#  from  -4* C  to  0°C.  is  where 
only  aggregation  is  observed.  Except  for  spiral  2#  the  ASD 
data  did  not  reach  such  a  high  temperature  and  both  N,  and 
X  decrease  between  the  last  two  loops  (18  and  19.  between 
temperature  -2*  C  and  0*0. 

It  is  revealing  that  all  observed  values  indicate  that 
the  slope#  A  .  always  terminates  at  around  10  cm'1  .  This  is 
in  accord  with  the  ASD  results  which  suggest  collisional 
breakup.  No  matter  what  the  height  and  environmental 

conditions  are#  the  observed  N0 values  are  within  three  orders 
of  magnitude  and  the  A  values  are  within  one  order  of 
magnitude. 

2.8  Counter-examples 

The  results  from  four  spirals  flown  in  1981  (on  an  AFGL  C-130 
aircraft)  are  shown  in  Figs.  2.29,  2.30,  2.31  and  2.32,  two  each  on  April 
14,  1981  and  on  April  17,  1981.  The  results  for  the  two  spirals  on  April 
14  do  not  indicate  the  three-stage  snow  growth  processes.  The  flights 
took  place  in  pre-frontal  showery  conditions  and  the  temperature  for  these 
spirals  is  rather  high.  The  changes  of  N^-X  for  the  four  cardinal  quadrant^ 
of  the  two  spirals  indicate  that  there  is  no  horizontal 
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Fig. 2. 30:  N  -A  trajectory  of  second  spiral  on  14  April,  1981 
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X  trajectory  of  first  spiral  on  17  April,  1981 


homogeneity.  Fig.  2.  33  is  the  example  from  the  first  spiral 
on  this  day.  The  changes  of  Ne-)\  are  disorganized  and 
incoherent  for  the  four  quadrants.  In  conclusion!  it  can  be 
said  that  the  ASD  will  work  only  when  the  prec ip itation  is 
widespread  and  horizontally  homogeneous.  It  will  not  work  in 
showery  conditions. 


The  changes  for  N#  and  from  the  spirals  on  April  17# 
1981  (Figs.  2.  31  and  2.32)  have  one  peculiar  behavior  which  is 
not  found  in  any  other  spirals#  namely  that  the  slope#  A  > 
kept  on  decreasing  even  beyond  the  10  cm*1  until  melting.  The 
synoptic  conditions  where  the  flights  took  place  are  a  weak 
system  between  a  surface  ridge  and  a  surface  trough  with  warm 
air  advection  from  the  gulf.  From  the  temperature  and  dew 
point  taken  on  the  plane#  the  environmental  conditions  are 
clearly  subsaturated  (Fig.  2.  34).  The  changes  of  total 
particle  concentration  from  the  1-D  precipitation  probe  with 
height  for  the  two  spirals  are  shown  in  Figs.  2.  35  and  2.36. 
Except  for  the  first  few  loops  of  the  spirals#  the 
concentration  is  decreasing  (compare  the  change  of 
concentration  with  height  for  the  1980  spirals  in  Fig.  2.  37>. 
Also  plotted  on  the  figure  is  the  second  moment.  The  changes 
of  the  sum  of  the  second  moment  are  comparable  to  those  of 
the  total  concentration.  Since  the  mass  of  a  snowflake  is 
approximately  proportional  to  the  square  of  the  diameter#  the 
sum  of  the  second  moment  is  proportional  to  the  total  mass. 


trajectories  for  the  four  cardinal  quadrant 
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So.  a  decrease  in  both  the  total  concentration  and  the  total 
mass  implies  that  particles  are  evaporating  due  to 
subsaturation  (compare  Figs.  2.  37  and  2.38  for  the  1980  ASDs). 
In  contrast  to  the  1980  2-D  data-  the  2-D  images  for  these 
two  cases  indicate  an  absence  of  needles  in  the  third  stage. 
This  also  suggests  that  the  atmosphere  is  subsaturated. 

A  passible  explanation  is  that  the  number  of  large 
particles  is  being  maintained  bg  aggregation.  while 
collisional  breakup  generates  minute  particles  that  are  not 
detectable  bg  the  1-D  precip  probe  (Fig. 2. 39).  Under 
conditions  of  deposition,  these  small  particles  will  grow 
into  detectable  sizes.  but.  under  conditions  of 
subsaturation,  these  small  particles  not  onlg  cannot  grow 
into  detectable  sizes,  but  rather  theg  are  evaporated.  The 
200-Y  probe  can  onlg  detect  particles  in  a  limited  range.  As 
depicted  in  Fig.  2.  39.  it  could  be  that  the  size  distribution 
detected  bg  the  probe  decreases  in  both  the  'intercept'  and 
the  'slope'  and  so  the  slope  will  decrease  below  the  10 
cm’1  limit  observed  in  the  other  spirals.  Figs.  2.  40  and  2.41 
show  the  evolution  of  the  last  few  spectra  in  the  two 
spirals.  The  behaviour  of  the  two  spirals  which  are  flown  on 
the  same  dag  and  in  similar  environmental  conditions  are 
consistent.  If  this  hgpothesis  is  true,  then  it  can  also  be 
postulated  that  the  collisional  breakup  fragments  are 
composed  mainlg  of  small  particles. 


.39:  Possible  size  distributions  beyond  the  range  detectable  by  the  200 
probe  at  three  different  times 
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Fig. 2. 40:  Snow  size  spectra  of  the  last  few  loops  for  the  first  spiral  on 
1-7  April,  1981 


Fig. 2. 41:  Snow  size  spectra  of  the  last  few  loops  for  the  second  spiral  on 
17  April,  1981 
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CHAPTER  3 


THEORETICAL  ANALYSES 


The  ASO  data  presented  in  the  last  chapter  reveal  that 
snow  growth  goes  through  different  distinct  stages  and  the 
intuitive  explanation  suggests  that  the  stages  are 
deposition.  aggregation  and  collisional  breakup.  If  this 
conjecture  can  be  supported  bg  a  rigorous  theoretical 
treatment.  then  the  intuition  is  confirmed  to  be  valid.  In 
order  to  accomplish  this.  various  theoretical  models  have 
been  developed  and  tested.  In  this  chapter,  such  models 
incorporating  the  physical  processes  of  vapor  deposition, 
sedimentation  due  to  differential  particle  fallspeed. 
aggregation,  collisional  breakup  and  rime-splintering  are 
presented. 

The  first  type  of  model  has  a  two-parameter  particle 
size  distribution  (2.1).  The  physical  processes  represented 
are  deposition,  aggregation  and  collisional  breakup.  Two 
forms  of  collisional  breakup  are  used.  The  first  one  assumes 
that  the  fragments  generated  from  collisions  are  distributed 
in  a  negative-exponential  distribution.  The  second  one 


assumes  that  the  number  of  fragments  generated  from  each 
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collision  event  ie  a  constant,  regardless  of  the  size  of  the 
colliding  particles.  Computations  are  performed  using  the 
values  of  N#  and  "X  from  the  first  loop  of  the  1980  spirals  as 
initial  conditions.  The  results  from  the  models  compare  well 
with  the  observations  and  so  the  intuitive  explanations  of 
the  observations  are  well  supported  by  theory. 

The  second  type  of  model  has  a  three-parameter  particle 
size  distributon  (2.4).  Again,  two  forms  of  collisional 
breakup  are  used.  In  this  type  of  model>  the  results  from 
using  a  negative-exponential  distribution  of  fragment  sizes 
give  satisfactory  results  compared  to  the  observations.  But* 
the  results  from  using  the  other  fragment  size  distribution 
do  not.  From  the  results  of  both  the  two-parameter  models 
and  the  three-parameter  models*  it  is  found  that  collisional 
breakup  is  an  important  snow— growth  mechanism.  especially 
when  the  temperature  is  close  to  melting.  Another  result  is 
that  when  two  snow  particles  collide*  they  are  more  likely  to 
*99r*9*t*  than  to  break  up.  However*  when  they  do  break  up* 
they  generate  numerous  small  particles. 

The  effect  of  different  initial  particle  concentration 
on  the  evolution  of  size  spectra  is  studied  with 
two-parameter  models.  Increasing  the  concentration  without 
increasing  the  mass  causes  the  change  from  deposition  stage 
to  aggregation  stage  to  start  earlier  than  before*  but  the 
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change  from  aggregation  to  collisional  breakup  stages  is 
delaged. 

The  problem  of  rime-splintering  is  studied  next. 
Theoretical  models  with  the  processes  of  deposition/ 
aggregation  and  rime-splintering  are  developed.  The  results 
are  compared  with  the  observational  results  from  the  spirals. 

They  show  that  rime-splintering  is  not  likely  to  be  a 

/ 

predominant  mechanism  in  this  kind  of  precipitation. 

To  end  the  chapter*  the  effects  of  melting  on  the  change 
from  snow-size  spectra  to  rain  drop  spectra  are  discussed. 


3.  1  Basic  equations  and  assumptions 

If  f(x*h.t>dx  is  the  concentration  of  particles  with 
mass  between  x  and  x+dx  at  height  h  and  time  t»  then  the 
change  in  the  number  density  due  to  nucleation*  vapor 
deposition*  sedimentation  due  to  differential  particle 
fallspeed*  aggregation*  collisional  breakup  and 

rime-splintering  can  be  formulated  using  the  stochastic 
collisional  model  as  (Scott*  1968;  Drake*  1972;  Passarelli* 
1978b  and  Srivastava.  1978): 
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-*•  ~~l  /S<x  J  x  x"  >f  <x  '»  h.  t)f  <  x",  h.  t)K  ( x  x") 

C 1— q (  x  xM)Jdx'dx" 


-^f(x,h,t) 

~  -Jr<Cw-V(x,3f<x,h,t»  +  c(t,h,<(x>  -  -|f<xf(x,h.t» 
/  /* 

♦  2^f<*”X,'h't><!<lt,'h't>K<1,~x,'  x  '  ^  *“*  '*  x  '  >d*  ' 

-  f  <x,  h.  t  >/f<x  h,  t  )K(  x<  x  ' ) q<  x»  x')dx' 

■b  (3.  1 ) 

-*•  /s(xl«/.x“)f(x'.h.t)f(x,,,h.t)K(x',x") 

Cl-q(x,,x,,)]dx,dx" 

-  f(x,h,  t>/f(x',  h.  t)K<x,  x')Cl-q<x,  x')]d«' 

+  ^->  CV(x)-uln'f(x#  h<  t)dx 

O 

where  w  *  updraft  velocity 

V<  x )  «  terminal  fallspaad  of  particle  with  mass  x 

u  *  terminal  fallspeed  of  supercooled  droplets 

x  *  rate  of  change  of  mass  due  to  vapor  deposition 

^(x)  «  Dirac  delta  function 

c(t<h)  *  nucleation  function 

a  collisional  kernel  between  particles  of 
mass  x  '  and  x " 

q(x',x">  *  probability  of  aggregation  when  two 
particles  of  mass  x'  and  xn  collide 

S<x!x'»x")  *  number  of  fragments  with  mass  between  x 
and  x+dx  when  particles  with  mass  x'  and 
x"  collide  and  break  up 

D  ■  diameter  of  snow  particle  with  mass  x 

d  ■  diameter  of  supercooled  liquid  droplets 


n'  ■  number  of  supercooled  liquid  droplets 

E  a  collisional  efficiency  between  snow  particles  and 
supercooled  liquid  droplets 


The  first  term  on  the  right  hand 


side 


i  s 


the 


sedimentation  due  to  differential  particle  fallspeed.  The 
second  term  is  the  increase  in  concentration  due  to 
nucleation.  The  Dirac  delta  function  is  used  because 
nucleation  can  only  increase  the  concentration  but  cannot 
affect  the  mass  or  the  reflectivity  factor.  The  third  term 
is  the  change  of  particle  spectrum  due  to  vapor  deposition  on 
snow  crystals.  The  fourth  term  is  the  production  of 
particles  due  to  crystal-crystal  aggregation.  The  fifth  term 
is  the  depletion  due  to  aggregation.  The  sixth  term  is  the 
production  of  particles  due  to  collisional  breakup.  The 
seventh  term  is  the  depletion  of  particles  due  to  collisional 
breakup.  The  eighth  term  is  the  increase  in  concentration 
due  to  rime-splintering.  Again*  this  term  can  only  affect 
the  concentration  but  not  the  mass  or  the  reflectivity 
factor. 


The  fragment  size  distribution  is  constrained  by  the 
following  relationship: 


S( x ! x  ' *  x">  *  0. 


if  x>x'  +  x"  (3.  2> 


Physical 

fragment 

original 


ly>  this  means  that  wher.  two 
can  be  larger  than  the 
particles.  Also* 


/ 


xS(  x  I  x  «")dx. 


particles  collide*  no 
sum  of  the  mass  of  the 


x  '+x" 


(3.  3) 
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This  means  that  the  total  mass  of  all  fragments  must  equal 
the  total  mass  of  the  parent  particles. 


From  these  basic  equations*  the  moment  conservation 
equations  can  be  derived.  Denoting  the  total  particle 


concentration*  n(h.t).  by 


n(h*  t) 


J  f (x*  h. 


t  >  d  x 


(3.  4) 


and  the  vertical  flux  of  particle  concentration*  n^(h.t)*  by 


n^  ( h  *  t )  *  J 


V<  x  )f  <  x,  h.  t ) d x . 


(3.  5) 


the  zeroth  moment  conservation  equation  can  then  be  written 

as 

JjLiUl  »  2h *<££ 

3*  u. 

-  J  f  <x',  h,  t>f(x".  h,  t)K<x'.  xM>q<x',  x">dx 'dx" 

J0J»  (3.6) 

+  -j*  S( x  «  x  x"  )f(x',h.t)f(x“  *h,t)K(x'.  x") 

C  i-q(  x  ' »  x“)ldxdx'dx" 

- Jj  f<x'«h*t)f(x"*h*t)K(x'*  x "  >  C  1-q (  x  '»  x">3dx'dx" 

Denoting  the  total  particle  mass  by  X  and  the  vertical 
mass  flux  by 


W‘ 


xf  (  x<  h*  t ) d x 


V(x)xf(x*  h*  t )  d  x 


(3.  7) 


<3.  8) 


The  first  moment  conservation  equation  can  be  written  as 
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=  Jjj, (-  +  £  if  <x.  h,  t)di 


(3.  9) 


Denoting  the  second  moment  of  the  mass.  which  mag  be 
interpreted  as  the  radar  reflectivity  factor,  as  Z(h.t)  and 
the  vertical  flux  of  the  second  moment  as  Zj(h.t). 


xl  f  (x,  h,  t ) d x 


Zj  *  Jvu)tx  f  <x,  h,  t )  d  x 


(3.  10> 


(3.  11> 


the  second  moment  conservation  can  be  written  as 


-  if  ».*>«< 


M.  J, - 

-*■  /  /  X  'x"f  (x  h.  t>f  <x".  h.  t>K(x  x")q<x  x")dx  'dx“ 

-4  #n  -f 


(3.  12) 


+  ±fl  f  -*-S(xix  '.  xM)f  <x  h,  t>f  (x",  h.  t)K(x  x“> 

*  Cl-q(x  x">3dxdx  'dx“ 

-  J  f  x1  f  (x.  h,  t)f  (x  '.  h,  t)K<x,  x  ')Cl-q(x.  x  ')3dxdx  ' 


Since  snow  is  falling  downward  towards  the  ground,  it  is 
easier  both  analytically  and  conceptually  to  take  a  reference 
level  of  height  and  designate  it  as  h*0  (e. g. .  at  the  top  of 
a  snowshaft  or  at  the  beginning  of  an  ASD).  Then  h  increases 
downwards. 


Since  the  observations  are  done  in  a  widespread.  slow 
updraft  wintertime  snowstorm.  steady  precipitation  can  be 
assumed.  When  this  assumption  is  applied  to  the  theory.  all 
time  derivatives  and  dependencies  in  the  moment  conservation 
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aquations  are  eliminated.  Moreover.  the  convergence  of 
upward  flux  of  moisture  is  balanced  by  the  increase  in  the 
downward  flux  of  precipitation* 

m  JjaJaL-  (3.  13) 

JA.  )i l 

where  ^  -  saturation  vapor  density. 


Integrating  with  respect  to  height  gives.  if  the 
downward  flux  of  prec ip i tation  is  taken  to  be  positive. 

^  fv'*  "  ^  fa  i0  (3.14) 


where  the  subscript  o  indicates  the  initial  condition.  By 
the  equation  of  state  and  the  Clausius-Clapeyron  equation, 
the  saturation  vapor  density  as  a  function  of  height  is 


id 

lz 


I 

z 


>  3exp  ( 


LiJL 

RV 


(3.  15) 


where  L  *  latent  heat  of  sublimation 
R  *  gas  constant  for  water 
i  *  atmospheric  lapse  rate 
Tq  a  initial  temperature 
At  h  *  0. 


u 

n 


exp 


(3.  16) 


and  denoting 
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A  =  <3.  17) 

UTo 

yields. 

h  -  X/.  -  «  (iu<**L  -  *> 

For  ,  Ah  »  1 

V  «  «,'L  o-  18> 

So.  for  a  known  atmospheric  lapse  rate.  the  precipitation 
mass  flux  is  determined  by  the  updraft  velocity. 

Other  assumptions  that  are  made  in  order  to  simplify  the 
computations  include:  The  rate  of  growth  of  a  particle  due 
to  vapor  deposition. 

x  -  g  <h )D  (3.  19) 

where  g(h)  is  a  function  depending  on  height  and  temperature. 
The  mass-diameter  and  mass-terminal  fallspeed  relationships 
are 


(3.  20) 

aDl 

(3.  21) 

where  Of.  p  .  a  and  b  are  constants  depending  on  the  snowfall 
type.  The  collision  kernel  is  formulated  according  to  the 
geometric  kernel. 

K <  x  ' ,  x "  )  -  -^-<D'+D"  >*  EIV(D"  >-V(D"  )  l 
where  E  *  collision  efficiency. 


(3.  22) 


From  the  clastic  properties  of  ice  (Hobbs.  1974),  it  can 
be  assumed  that  the  larger  a  snow  crystal,  the  more  easily  it 
breaks  upon  collision.  Mathematically,  the  coalescence 
probability  can  be  represented  by 


q  (  x  ' »  x  "  > 


e 


(3.  23) 


where  c  is  a  constant.  The  exponential  form  is  chosen  to 
facilitate  the  mathematical  analysis.  A  plot  of  the 
coalescence  probability  is  shown  in  Fig.  3.1. 


All  the  mathematical  analyses  that  follow  make  use  of 
the  above  assumptions. 


3.  2  Two-parameter  formulations 

The  snow  size  spectrum,  f(x,h)dx,  can  be  represented  by 
a  Ounn  and  Marshall  (1958)  type  negative  exponential 
distribution, 

f  (x,  h)dx  -  N0  (h)e’A^JJ>  dD  (3.24) 

where  Ne(h)  and  ?\<h)  are  functions  of  height.  In  this 
formulation,  there  is  one  independent  variable,  h,  and  two 
dependent  variables,  N,  and  }\  .  (Kessler,  1969i  Srivastava, 


1971,  and  P 


Hi,  1978a,  b) 
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Fig. 3.1:  Coalescence  probability,  q,  as  a  function  of  the  total  diameter 
D'-H)"  and  the  constant  c 


Dlomefer  of  porficfo  >r>  cm 

Number  of  fragments  resulting  from  collision  of  a  2 
and  particles  of  various  sizes  for  diffemt  A  s 


mm  particle 
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3.2.1  Exponential  fragment  size  distribution 

Previous  studies  (e. g. •  Var diman.  1978)  have  shown  that 
the  larger  the  crystals.  the  more  fragments  are  generated 
during  collisions.  The  laboratory  investigations  which  are 
to  be  described  in  appendix  2  indicate  that  the  fragments  are 
in  an  approx imately  negative  exponential  size  distribution. 
These  motivate  the  formulation  of  the  fragment  size 
distribution  to  be 

S (  x  !  x  ' <  x  “ )  -  (x  '+x">  Ale'A*  (3.25) 

where  A  is  a  constant.  The  total  number  of  fragments  is 

js<  x  1  x  '»  x“)dx  -  (x'  +  x“>/\  (3.26) 

which  implies  that  the  larger  the  particles.  the  more 
numerous  are  the  fragments.  Plots  of  the  number  of  fragments 
resulting  from  the  collision  of  a  2  mm  particle  and  particles 
of  various  sizes  for  different  /\'s  are  shown  in  Fig.  3.  2. 


dXdY. 


*  fiLft-JLAf' - 1  /  xf  Y*  (X+Y)1  JX^-Y1  !e  * 

t>+ ob^*4  4-c 

Then  the  second  moment  conservation  equation  is 

Hi-  2*  gN.  Pc  (3  +2)  2  «*  -j-Eal, 

iJL  r  OUc)^'V 


♦  4  H  N,  -y~Ea !». 

A  V*?*4 

-  2rf*N*  j-Eal, 
yVfT^t 


4x  N0L  -j-Ealt 


A  l\+4W*  + 

2  j-Eal, 

a*o^*4 


(3.  28) 


where 


I(  =  fjfx^Y^  (X+Y)1"  <Xl-Yfc  )e'f<+'*>  dYdX 
Ix  >  J-j <x+Y)*  IX^-Y^I#'**^  dXdY 
I3  «  ~jj  K2?  (X+yf  !Xfc-Y kle'^>  dXdY 


dXd  Y 


(3.  29) 


(3.  30) 


<3.  31) 


and  where  P  is  the  gamma  function.  The  integral  can 
expressed  as  hypergeometric  functions#  e.  g.  . 

I  *  r<b+2  +4) 

F  <l,b+2f+4f  p+Zi-lfi  XF  <1.  b+2j?+4i  b+  /»  +2;  --) 


(3.  32) 

F  (1,  b+2j>+4/  p+3 i-J-)  7lF  <1.  b+24+4 ib+6+3i-}-) 


F;  <  1#  b+2/J+4» 


F(  <1.  b+2p+4>  b+/?+4 


Yi 


Making  usa  of  the  vertical  flux  of  precipitation  mass. 


the  dependant  variable  N  can  be  removed  via 


N 


o 


a.ot  r(b*p+o 


(3.  33) 


The  function  for  the  rate  of  change  of  particle  mass  due  to 
vapor  deposition  can  be  avoided  by  making  use  of  the  first 
moment < 


g(h) 


(J rci+H) 


(3.  34) 


Finally<  the  change  of  with  height  is  given  by  an 
implicit  equation 

5L.  „  ffl>  , 

SAT  ~Xf  f  0 

af e  X(  », _ 

a b< p (\+c )**a/,t*  r<b+pn  r<b+2^-n) 


A  a*'jJ  P(b+{*  +1)  r<b+2^+l) 

4<4>ex,  ,t 

/\a*l(l  (  X+c)  **P+4 

f  2fE*»  ^  h 

a  4^  f(b+  (l  +1 )  f<b+2^  +  l> 


(3.  39) 


f<b+fl>  P<b+2p 
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_  xiSh  tan  h _ 

a*  p  <  P( b+p+1)  r<b+2^+l) 

\  is  obtained  by  numerically  solving  (3.35)  in  height  steps. 
Once  \  is  obtained*  N.  can  be  found  from  (3.33). 

# 

In  the  computations*  the  values  of  N(  and  )i  from  the 
first  loop  of  an  ASD  are  taken  as  the  initial  conditions. 
Because  of  the  lack  of  published  values  on  the  collision 
efficiency*  E*  it  is  taken  to  be  1.0  for  all  cases.  The 
magnitude  of  the  constants  a*  b*  <X  and  are  guided  by  the 
values  of  Locatelli  and  Hobbs*  1974.  But  the  actual  values 
of  a*  b*  C*  «  |3  *  c  and  A  used  in  the  computations  are 

determined  by  trial  and  error  so  that  the  theoretically 

predicted  values  duplicate  the  values  from  observations. 
This  approach  for  determining  the  constants  may  seem  rather 
empirical  but  due  to  the  lack  of  published  values  of  these 

constants*  this  is  the  only  way.  In  fact*  one  of  the 

purposes  of  this  study  is  to  compare  the  theoretical  values 
with  the  observational  findings  so  as  to  deduce  the  relative 
magnitudes  of  the  various  constants. 

The  results  from  this  formulation  are  plotted  against 
the  observations  in  N. -A  space.  Figs.  3.  3  and  3.4  show  the 
results  for  25  Feb.  *  1980  and  26  Feb. *  1980.  The  values  of 
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.wo-parameter  model  evaluations  tor  spiral  on  Zb  February, 
>reakup  fragment  size  distribution.  See  text  for  details. 


26-80 


cm/s*  so  in  has  been  sat  to  10  cm/s  in  all  the  calculations. 
The  effect  of  different  values  of  w  on  the  computation  is 
being  shown  later.  The  results  from  the  theory  also  show 
that  snow  growth  goes  through  three  distinct  stages.  In  the 
first  stage*  the  intercept*  N0*  increases  rapidly  while  the 
slope*  *)\  <  remains  constant.  This  dupicates  the  deposition 
stage  in  the  observational  data.  The  second  stage  of  growth 
is  where  both  N(  and  )\  decrease  rapidly.  This  duplicates  the 
aggregation  stage  from  the  observations. 

The  results  from  the  theoretical  model  but  with  c  set  to 
0  (i.e.  no  collisional  breakup  process)  are  also  plotted  in 
Fig.  3.  3  as  a  dashed  line.  The  results  from  this  formulation 
have  only  two  stages.  The  first  stage  is  when  increases 
while  ^  remains  almost  constant.  The  second  stage  is  when 
both  N,  and  ^  decrease  rapdily.  In  this  formulation*  there 
seems  to  be  no  end  to  the  decrease  of  and  7\  . 

The  results  from  the  theoretical  model  with  collisional 
breakup  (c>0>  have  a  third  stage  where  N  increases  and 

w 

^remains  constant.  From  previous  studies  (e.  g.  »  Srivastava* 
1978)  it  is  concluded  that  aggregation  and  collisional 
breakup  processes  can  balance  one  another  so  as  to  produce  an 


Fig. 

Date 

a 

b 

a 

8 

E 

V 

c 

A  Ns 

3 

2-25-80 

100 

0.15 

0.002 

2.15 

1 

10 

0.055 
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4 

2-26-80 

100 

0.15 

0.002 

1.55 

1 

10 

1.0 

3500 

7 

3-08-80 

100 

0.15 

0.002 

0.75 

1 

10 

0.65 

10000 

8 

3-08-80 

100 

0.15 

0.002 

2.0 

- 

10 

1.5 

9 

2-25-80 

100 

0.15 

0.002 

1.95 

1 

10 

0.07 

1750 

10 

2-26-80 

100 

0.15 

0.002 

1.25 

1 

10 

0.25 

15 

11 

3-08-80 

100 

0.15 

0.002 

0.75 

1 

10 

0.65 

20 

12 

3-08-80 

100 

0.15 

0.002 

2.0 

1 

10 

1.5 

0. 1/0.5 

Table  3.1:  Values  of  the  physical  parameters  used  in  the 
two-parameter  model  (c.g.s.  units) 


a 

s 

b 

c 

A 

X 

.  .  _o_. 

X  in 
model 

0.002 

2.15 

0.15 

0.055 

9000 

9.99 

8.30 

0.002 

1.5 

0.15 

1.0 

3500 

16.24 

11.35 

0.002 

0.75 

0.15 

0.65 

10000 

9.99 

7.43 

0.002 

2.0 

0.15 

1.5 

1750 

-0.75 

8.11 

Table  3.2:  Tabulation  of  X 
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equilibrium  slope.  This  is  clearly  the  case  here.  The 
introduction  of  the  collisional  breakup  process  into  the 
model  halts  the  decrease  in  X  and  leads  to  an  equilibrium. 
The  reason  that  N.  keeps  on  increasing  is  because  vapor 
deposition  is  still  going  on  and  if  X  is  constant*  N0 has  to 
increase  so  as  to  accomodate  the  increase  in  mass. 

Fig.  3.  3  shows  the  effect  of  varying  the  updraft  velocity 
while  keeping  all  other  parameters  constant  for  the  case  of 
26  Feb.  »  1980.  The  differences  in  updraft  velocities  affect 
only  the  position  where  the  first  stage  changes  to  the  second 
stage  and  does  not  affect  the  final  X-  Since  the  updraft 
velocity  affects  the  deposition  rate.  so  a  change  of  updraft 
velocity  will  naturally  affect  where  the  first  stage  changes 
to  the  second  stage.  The  stronger  the  updraft,  the  greater 

N  has  to  be  in  order  to  accomodate  the  additional  mass.  The 

o 

fact  that  the  updraft  does  not  affect  the  final- X  is 
consistent  with  former  studies  of  Srivastava  (1978)  that  the 
equilibrium  X  resulting  from  aggregation  and  collisional 
breakup  is  independent  of  updraft  velocity.  The  observations 
from  section  2.8  show  that  for  subsaturation  cases,  the  final 
slope  will  be  affected.  This  might  be  due  to  the 
microphysical  processes,  as  are  explained  in  that  section. 

The  results  from  the  case  of  8  (larch.  1980  are  also 
plotted  in  Figs.  3.  6  and  3.7.  This  is  the  case  where  there  is 


26-80 


[ill] 


not  an  initial  stag*  where  N,  increases  while  a  remains 
constant.  In  the  process  of  finding  the  set  of  parameters  in 
which  the  theoretical  model  can  best  duplicate  the 
observations*  it  is  discovered  that  two  sets  of  parameters 
can  more  or  less  duplicate  good  fits  (Table  3.1).  One  of  the 
sets  requires  |3  .  the  mass-diameter  power  to  be  2  which  is 


within  published  values.  The  other  set  requires  p  to  be  0.79 
which  is  too  small.  Moreover*  this  set  requires  A  *  the 
breakup  exponent  to  be  10000*  which  implies  the  number  of 
breakup  fragments  to  be  high.  Results  from  the  other  breakup 
formulation  and  from  the  three-parameter  formulations  can  be 
used  to  clarifg  the  situation. 


A  simple  analgtical  expression  of  the  equilibrium  slope 
can  be  obtained  bg  assuming  no  vapor  deposition.  Then  the 
first  moment  conservation  equation  becomes* 


-\-£-  - 0 

which  implies 


V  '  ""***"* 


(3.  36) 


The  second  moment  equation  is 


_df  -  N)of-|-EaI, 
>h  (  ^c)^ 


Ne2«(^-Eall 

aV7^ 


A  <  A  +c  I1**'* 


<3.  37) 


>rv  <  x  +c  r2?** 


From  these  two  equations  and  setting 
equilibrium  slope.  o.  is  expressed  by 


to 


O. 


the 


n,\!.  i,  c  -j,  *  x. 

A  a  J  • 


2  Ix_ 


i 


(3.  38) 


This  expression  shows  that  the  equilibrium  slope  is  dependent 
on  the  density  (  c*  and  ji  ).  the  fallspeed  (b>  and  the  breakup 
parameters  (  /\and  c)  and  is  independent  of  the  mass  flux,  the 

updraft  and  the  collisional  efficiency.  X  obtained  by 

o 

using  the  values  of  parameters  used  in  the  model  agrees  well 
with  the  final  A  value  obtained  in  the  model  (Table  3.2). 
The  only  exception  is  the  8  March.  1980  case  with  in 
which  the  obtained  is  even  unrealistically  negative.  As 
will  be  seen  in  the  later  sections,  computations  using  mZ 
for  the  8  March.  1980  case  fail  to  produce  satisfactory 
results,  which  says  something  about  the  suitablity  of  this 
set  of  parameters.  (3.38)  indicates  that  the  equilibrium 
slope  depends  only  on  the  density  of  snow.  the  coalescence 


probability 

and  the  breakup 

fragment  distribution. 

The 

observations 

show  that 

the 

equilibrium 

slopes 

are 

approximately 

the  same 

for 

all  flights. 

So.  if 

the 

equilibrium  slopes  are  approx imately  the  same  for  the  various 
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flights<  this  might  suggest  that  the  density  of  snow  and  its 
breakup  properties  are  interrelated. 

The  values  of  N  and  where  the  first  stage  of  snow 

c 

growth  changes  to  the  second  stage  can  be  analytically 
obtained  by  neglecting  collisional  breakup.  Then  by  setting 

-&■ to  6*  0< 

{JL£dLlLkti&}  £jf£r 

h 


c  gitxp  ro^fo  { 

u+p+o  t  £  I,  h 


<3.  39) 


.  Cl 


pf 


The  values  of  N,  and  }i  obtained  from  (3.  39>  for  the  different 
spirals  are  shown  in  Table  3.3.  These  values  are  close  to 
the  values  obtained  from  the  model.  But  this  shows  that  the 


collisional  breakup  process  does  have  a  weak  effect  on  the 
change  from  the  first  stage  to  the  second  stage. 
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No 

1.23  0.4662 

0.18309 

0.075261 

X 

5.38  38.9 

28.5 

21.2 

Table  3.3:  First 

turning  point 

of  N0 

Fig. 

Date 

a 

b 

a 

0 

E 

w 

c 

A 

No 

13-15 

2-25-80 

100 

0.15 

0.002 

2.1 

1 

10 

0.1 

5000 

16-18 

3-08-80 

100 

0.15 

0.002 

0.75 

1 

10 

0.5 

10000 

19-21 

3-08-80 

100 

0.15 

0.002 

2.0 

1 

10 

3.0 

8000 

22-24 

2-25-80 

100 

0.15 

0.002 

1.7 

1 

10 

0.2/0. 1 

25/7.5 

25-27 

3-08-80 

100 

0.15 

0.002 

0.75 

1 

10 

0.5 

5 

Table  3.4:  Values  of  physical  parameters  used  in  the  three-parameter  model 

(r.g.s  units) 
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In  order  to  illustrate  the  effect  of  changing  the  values 
of  the  various  physical  parameters  on  the  results*  the 
results  from  computations  using  different  values  of  the 
parameters  are  plotted  in  Fig.  3.  8.  The  case  used  is  the 
spiral  on  25  Feb. *  1980.  The  values  of  c  used  ranges  from 
0.05  to  1.  5.  which  means  that  q<  the  probability  that  two 
particles  will  aggregate  when  they  collide*  is  above  80%  for 
most  particles.  The  values  of  A  used  range  from  1750  to 
10C00.  A  >5000  implies  that  the  number  of  breakup  fragments 
can  easily  be  more  than  10.  From  the  values  of  c  and  A  *  it 
appears  that  when  two  snow  particles  collide*  they  are  more 
likely  to  aggregate  than  to  break  up*  but  when  they  do  break 
up.  they  generate  numerous  fragments.  The  larger  the  parent 
particles,  the  more  numerous  the  fragments. 


3.2.2  Constant  fragment  number  distribution 


In  this  formulation,  the  fragment  size  distribution  is 
assumed  to  be 


S  (  x  S  x  ' «  x  “  >  = 


Nj  U,  -  >. 


(3.  40) 


This  means  that  when  two  snow  crystals  collide  and  break  up* 


they  generate  Ns equal  mass  fragments.  So*  the  number  of 
fragments  does  not  depend  on  the  size  of  the  parent 


25  -80 


distribution 
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particles.  but  the  mats  of  fragment  does. 


Using  a  technique  similar  to  that  of  the  previous 
section.  the  change  of  A.  the  spectrum  slope,  with  height 
can  be  written  as  an  implicit  equation. 


rcwijito 


a-j~E  If  I, _ 

N^a «( |)  f(b+^+l>  r<b+2|i+l) 


(3.  41) 


(-ujjj-  a-|e  x,  A**!*"  h 

aWff(b-»-  (1  -»-l>T<b+2f*  +1) 

i  (-**&•  M6*.  (i^v 

a<rfp(V* €)***?**  T<b+^  ♦!>  P(b+2^J +l> 


where  I,  and  I3  are  (3.29)  and  (3.31)  in  the  previous  section. 
Again.  the  implicit  equation  is  solved  numerically  in  height 
steps  and  NQ's  are  obtained  from  (3.33).  Results  from  this 
formulation  for  the  three  different  cases  are  plotted  against 
observations  in  Figs.  3.  9-3.  12  and  the  values  of  the  physical 
parameters  used  are  listed  in  Table  3. 1.  Except  for  the  case 
of  8  March.  1980.  the  results  from  this  formulation  are  very 
similar  to  the  results  from  the  exponential  fragment  size 
distribution  formulation.  Fro*  the  quality  of  fit  to  the  29 
Feb.*  1980  and  the  2A  ..  *80  cases,  it  appears  that  there 
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is  not  much  difference  between  the  two  f ormulations. 
However,  for  the  8  March  case,  the  model  from  this 
formulation  can  fit  the  observations  only  by  using  the  set  of 
parameters  which  require  the  mass-d iameter  power  to  be 

0.79.  When  /;>  equals  2,  the  model  results  cannot  duplicate 
the  observe*:.'.  oVis  even  after  many  trials  (e.  g.  ,  Fig.  3.  12). 
Another  point  is  the  fragment  number,  N* »  required  in  the 
model  with  (i  *0.  75  is  only  7.5  which  is  not  high.  A  full 
interpretation  of  the  implications  from  the  values  of  the 
physical  parameters  will  be  discussed  after  the  presentation 
of  the  three-parameter  models. 


3.  3  Three-parameter  formulations 


In  this  formulation,  snow-size  spectrum  is  represented 


by 


f  <  x , h  ) d x  »  N^(h>  e  *  D  dD 


(3.  42) 


where  N^,,  A*  and  <r  are  functions  of  height,  h.  In  this 
formulation,  h  is  the  independent  variable,  while  N»,  A*.  and 
rare  three  dependent  variables. 


i\t  il 


Jm 
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3.3.1  Exponential  fragment  size  distribution 


The  fragment  size  distribution  S(xix',x")  is  formulated 
according  to  <3.  25)  as  <  x  '-n"  )  AVa*  .  From  this  formulation, 
the  zeroth  moment  conservation  equation  can  be  written  as 


3"  N l  -$-EaIA  2N^  -y-Eal^ 


<  Vc  ** 


*AN*  -^"Ea  Is 

\**1*x*‘+l 


Ad  N*  -jUdl. 


where 


l)  «  I Y^tX+Y)1  <Xk-Yl  )e'(<*^  dYdX, 

!r  -  f  fxf**"  Yr<X+Y>1  lXk-Y4l  e'cKf^  dXdY 
The  first  moment  conservation  equation  is 

?JL  *  Xf* 

The  second  moment  conservation  equation  is 


(3.  43) 


(3.  44) 


(3.  45) 


(3.  46) 
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-  - - +  7^7. 


20UM^J-EaIr 
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j-e«V 

A(  X*c>k>,''“''“' 
2d*N* 

<  X,«  >k**4"* 


(3.  47) 


Ill 


where  ^ 

*  J  JX?+<r  (X+Y)1'  (Xt-Yt)e'°‘^>>  dYdX  (3.48) 

I1  ■  f"*f  X^Y*  (X+Y)4  SX^-Y1  ie'^  dXdY  (3.49) 

•  • 

Similar  to  that  of  the  two-parameter  formulation*  can 
be  removed  by  using  the  prec ip itation  mass  flux* 
v 

N  -  -&-£■*■ - r  (3.  50) 

*  a*  ra+p+w) 

and  by  making  use  of  the  first  moment>  the  rate  of  change  of 
particle  mass  due  to  vapor  deposition  can  be  avoided* 

r  (L+ptr+i) 

g  ( h  )  *  — ■“ -  (3.51) 

After  these  manipulations*  two  implicit  equations  in  X*  and 
(rare  derived  as  functions  of  height. 

rr k+T  ,  ,  2<T 

“j*  =»  - - - - - - —Cl  xF  lnx  e  d x ]--- 

(i  Xv  (3r(lKT+i>  J,  M. 

♦  ------  C [k^P*  In  x  e~*  dxl-~ 

+  ilh  u _ 
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*  Ua  ,  (3  52) 
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(3.  53 > 


The  integrals  /  X  lnXe  dX<  /  X  r  lnXe  dX  and 

I  X  (  lnxe'*  dX  art  computed  using  a  Gauss— Laguerre 
quadrature.  Bg  using  a  second-order  Runge-Kutta  algorithm. 
X,and  0"  are  computed  in  height  steps  from  this  set  of  two 
equations.  The  N^<  X^and  at  the  beginning  of  the  ASD  data 
are  taken  to  be  the  initial  conditions  for  computation.  The 
choices  of  a.  b<  c«  ft  <  »  &  and  f\  are  guided  bg  the  results 

from  the  two-parameter  formulations. 


Figs.  3.  13.  3.  14  and  3.  19  compare  the  model  results  of 
Nj.  X*  <T  against  the  observations  from  25  Feb..  1980.  In 


order  to  show  the  effects  of  changing  the  values  of  the 
phgsical  parameters.  results  from  the  various  computations 
are  also  plotted  on  the  same  diagrams.  For  N^.  the  model  can 
djelicate  the  three  stages  in  the  observations.  First,  there 
is  a  slow  increase,  then  a  rapid  rise  and  a  rapid  fall  and 
finallg  Ng takes  on  a  more  or  less  constant  value.  Just  like 
the  from  the  observation.  the  Xj.  from  the  computation 
decreases  monotonicallg  with  height.  and  then  becomes 
approKimatelg  constant.  For  (T  .  the  model  result  is  negative 
for  a  short  while.  Then  it  changes  rapidlg  to  become 
positive  and  then  becomes  constant.  Comparing  this  with  the 
observational  results.  the  change  from  positive  to  negative 
occurs  too  earlg  for  the  theorg.  The  positive  values 
attained  bg  the  model  are  much  too  high.  In  fact  it  appears 
that  the  model  result  overshoots  the  equilibrium  value  before 
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it  decreases  beck  to  the  final  value.  In  summary  the 
theoretical  results  can  duplicate  the  different  stages  of 
snow  growth.  However,  the  quantitative  agreement  between  the 
observations  and  the  theory  is  not  as  good  as  in  the 
two-parameter  formulations.  There  are  three  reasons  for 

this.  First,  a  three-parameter  analysis  of  data  is  more 
sensitive  to  noise  in  the  data  than  a  two-parameter  data 
analysis.  Secondly,  it  is  always  easier  to  fit  two  variables 
than  to  fit  three  variables.  Finally,  the  accuracy  of  the 
smallest  few  size  channels  of  the  200-y  are  questionable. 
Ideally.  if  the  data  from  the  cloud  probe  could  supplement 
this  gap.  the  comparisons  between  the  data  and  theory  would 
be  better. 

Figs. 3.  16-3.  21  show  the  results  from  this  formulation 
against  the  observations  on  8  March.  1980.  Guided  by  the 
results  from  the  two-parameter  model,  two  sets  of  parameter 
values  are  attempted  and  they  are  listed  in  Table  3.  4.  The 
N^and  from  the  set  of  parameters  are  similar.  Both  N^and 
X^decrease  rapidly  and  then  become  constant.  These  basically 
agree  with  the  observational  results.  In  the  observations 
for  N^.»  except  for  the  third  and  fifth  loops,  it  decreases 
monotonically.  The  theory  does  not  have  the  'bumps'  as  such. 
It  is  not  clear  that  the  increases  in  N^at  these  two  loops 
are  real  or  simply  noise.  The  final  value  attained  by  in 
the  theory  is  higher  than  that  in  the  observations.  The 


Na  from  three-parameter  model  evaluations  for  spiral  on  8  March,  1980,  exponential 
breakup  fragment  size  distribution,  6=0.75 
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major  difference  between  the  two  sets  of  parameters  lie  in 
the  results  for  .  The  results  from  the  set  that  requires 
p>  =0.  75  decreases  with  height  and  then  becomes  constant  while 
the  results  from  the  set  that  requires  |3*2.  O  decrease  more 
rapidly  and  then  increase  again  to  become  positive.  It 
appears  that  the  set  of  parameters  that  0.75  duplicates 
the  observations  better.  However#  for  ^ »  the  mass-d iameter 
power,  to  be  0.  75  is  a  bit  too  low.  This  implies  the  density  of 
snow  is  small.  The  2-D  images  do  not  indicate  anything 
particularly  different  between  this  case  and  the  cases  on  25 
Feb.  #  1980  and  26  Feb. <  1980.  The  values  of  the  breakup 
parameters  imply  that  the  number  of  fragments  is  high#  which 
is  expected  because  <T  changes  from  positive  to  negative 
meaning  small  particles  are  generated  faster  than  large 
particles. 


3.3.2  Constant  fragment  number  distribution 

In  this  formulation>  the  fragment  size  distribution  is 
assumed  to  be 

S  (  x  i  x  '  #  x  "  )  =  N.  ^<x - - >  (3.54) 

*  Ns 

Using  a  technique  similar  to  that  of  the  previous  section#  a 
set  of  two  implicit  equations  are  obtained# 
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where  l^i  I^and  I-^are  from  (3.44).  (3.48)  and  (3.49)  in  the 
previous  section. 

The  results  from  these  formulations  using  the  2S  Feb. • 
1980  case  as  initial  conditions  are  shown  Figs.  3.  22.  3.23  and 
3.24.  Even  after  numerous  attempts  at  adjusting  the 
parameter  values.  N+<  A^and  <r  obtained  from  this  formulation 
cannot  be  made  to  fit  the  observations.  As  a  last  resort, 
one  set  of  parameters  is  used  from  the  beginning  until  the 
maximum  of  N^and  then  a  second  set  of  parameters  is  used  from 
then  onwards  until  the  end  of  the  computations.  Even  after 
all  these  efforts.  N^«  and  <T  still  cannot  fit  the 
observational  results  desirably.  This  probably  implies  that 
the  exponential  fragment  size  distribution  is  a  better 
breakup  formulation  than  the  constant  fragment  number 
distribution.  Physically,  it  is  more  reasonable  to  assume 
that  when  particles  collide,  they  do  not  generate  a  constant 
nunmber  of  fragments.  but  rather  the  larger  the  parent 
particles.  the  more  numerous  the  fragments.  Based  on  the 
laboratory  experiments  described  in  appendix  2.  it  is  also 


more  reasonable  to  assume  that  the  particle  sizes  are  in  a 


negative  exponential  distribution  than  to  assume  that  they 
have  uniform  sizes. 

Figs.  3.  23/  3.26  and  3.27  show  the  results  from  this 
formulation  against  the  observations  from  8  March/  1980. 
Only  the  set  of  parameters  having  0.  73  is  being  used. 
(The  values  of  the  physical  parameters  used  are  listed  in 
Table  3.4.  )  Both  N^,and  duplicate  the  observations 
competently.  However/  although  CT  decreases  initially  from 
positive  to  negative/  it  fails  to  maintain  a  constant  value. 
Instead/  it  gradually  increases  to  become  positive  again. 
This  behavior  is  similar  to  the  result  from  the  exponential 
fragment  site  distribution  with  p  equaling  2.  Considering 
the  results  from  both  spirals/  it  appears  that  a  negative 
exponential  size  spectrum  is  the  better  model  for  the  breakup 
fragments. 


3.4  Effect  of  different  initial  concentrations 

In  order  to  study  the  effect  of  seeding  (either  natural 
or  artificial)/  an  attempt  is  made  to  modify  the  initial 
values  of  parameters  so  that  the  total  concentration 
increases  but  not  the  total  mass.  For  the  two-parameter 
formulation/  assuming  the  mass-d iameter  relationship  to  be 


(sj;un)  s'fi  o 


X 


C*Dr  *  the  total  mass  is 


Jo (D^N,  e'^ 


Constant  total  mass  means 


N  =  constant 
0 


<3.  57) 


(3.  58) 


So.  the  change  of  concentration  without  the  change  of  mass  is 


where  N0l 


(3.  59) 

(3.  60) 


Using  the  case  of  26  Feb.  »  1980  as  an  example.  the 
initial  total  concentration  is  increased  5  and  10  times, 
without  changing  the  total  mass.  The  computations  are  done 
using  the  same  set  of  parameters  as  is  done  before  in  section 
3.2.2  and  the  results  are  plotted  in  Fig.  3.  28.  The  results 
indicate  that  the  change  in  the  initial  total  concentration 
affects  the  change  from  the  deposition  stage  to  the 
aggregation  stage  but  does  not  affect  the  final  equilibrium. 
Increasing  the  initial  concentration  5  times  causes  the 
aggregation  stage  to  start  approx imate ly  50  meters  earlier 
than  before  but  delags  the  collisional  stage  bg  approx imatelg 
500  meters.  Increasing  the  initial  concentration  bg  10  times 
causes  the  aggregation  stage  to  start  approx imatelg  100 
meters  earlier  and  delags  the  collisional  stage  bg  600 


meters. 


cm 


Physically/  the  increase  in  concentration  without 
increasing  the  mass  means  that  the  number  of  collision 
increases.  Quantitatively/  the  number  of  collisions  is 
represented  by 

j  (X+Y)1-  !Xt-Yl'i  e~(t*V  d Xd Y 

Since  Nv  's  and  ^'s  are  related  by  (3.60)/  higher  values  of 
N( and  lower  values  of  'X  means  higher  number  of  collisions/ 
resulting  in  an  increase  in  aggregation  events.  Therefore/ 
the  aggregation  stage  starts  earlier.  On  the  other  hand/ 
increasing  the  concentration  without  increasing  the  mass 
means  that  the  average  size  of  particles  are  smaller.  So/  it 
would  take  longer  for  them  to  grow  to  such  a  size  that 
collisional  breakup  becomes  efficient.  Since  neither 
aggregation  nor  collisional  breakup  affects  the  total  mass/ 
so  the  final  equilibrium  position  between  these  two  processes 
depends  only  on  the  slope/  'A  z  and  not  on  the  initial 
concentrations. 


3.9  Investigations  of  rime-splintering 

Mossop  (1976)  concludes  from  laboratory  experiments  that 
secondary  ice  particle  production  due  to  rime-splintering 
depends  on  the  number  density  of  supercooled  drop  greater 
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than  24  microns.  Approximately)  one  ice  splinter  is  thrown 
off  for  the  accretion  of  every  250  drops  greater  than  24 
microns.  The  temperature  range  for  occurrence  is  between 
-3*0  and  -8* C  with  the  production  rate  being  greatest  at 
-5*C.  Fig.  3.  29  shows  the  number  density  of  particles  greater 
than  24  microns  as  measured  by  the  Axially  Scattering 
Spectrometer  Probe.  Based  an  Fig.  3.  29  the  splinter 
production  rates  of  the  three  spirals  can  be  estimated  by 
assuming  that  these  are  supercooled  water  drops. 

In  unit  time*  the  number  of  accreted  drops  by  a  single 
snow  particle  is 

7[  E(~*  ■»■  -^■>1CV<D)  -  V(d ) 3n '  (3.61) 

where  D  =  diameter  of  ice  particle 

d  *  diameter  of  supercooled  droplet 

V(D>.  V(d)  *  terminal  fallpeed  of  snow  particle  and 
drop  respectively 

n'  *  number  density  of  supercooled  droplets 
E  =  collisional  efficiency 

Representing  the  number  density  of  snow  particles  with 
diameters  between  D  and  D+dD  as  N(D>dD>  then  the  total  number 
of  accretion  events  per  unit  time  per  unit  volume  is 

Ju  E(-j  +  CV(D)  -  V<  d  >  3n 'N(D) dD  (3  6 7>) 

dD,  V  *  aDV,  and  V(d)«V(D),  then  the 


Assuming  N(D>dD  *  No  e 
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integral  becomes 

r  ru+3)  .  rci>+2)  dL  ra>*o  ,  _ 

I  "  -7^r '  >>“  - 

Using  typical  values  of  a  and  b  (Locatelli  and  Hobbs.  1974> 
and  the  observed  values  of  n'  from  Fig.  3.  29.  the  accretion 
rates  for  each  spiral  can  be  estimated.  Based  on  Mossop's 
conclusion  that  one  splinter  is  produced  for  every  250 
accretions.  the  splinter  production  rate  can  also  be 
estimated.  The  accretion  rates  and  the  production  rates  for 
the  three  spirals  are  listed  in  Table  3.5. 


The  particle  depletion  rate  due  to  aggregation  can  be 
estimated  using  the  stochastic  model  with  a  geometric  kernel. 


iP 


N(D  )N(D 

“I  O  Z 


>-$*• 


!D,‘-d‘ 


(D, 


+v 


dD  dD„ 


(3.  64) 


The  estimated  depletion  rates  for  the  three  spirals  are  also 
listed  in  Table  3.  9. 


From  Table  3.5  it  appears  that  the  splinter  production 
rate  is  compatible  with  the  particle  depletion  rate.  If  the 
particles  detected  by  the  ASSP  are  really  supercooled 
droplets.  then  the  Mossop-Hal let  splintering  mechanism  could 
account  for  the  secondary  particle  production.  As  mentioned 
before.  the  ASSP  probe  operates  on  a  light  scattering 
principle  and  can  easily  give  false  counts  from  the  light 
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scattered  or  reflected  of  the  snow  particles.  Comparing 
Figs.  2.  9.  2.11  and  2.13  with  Fig.  2.  38.  the  ASSP  concentration 
and  the  1-D  precipitation  probe  concentration  increase  and 
decrease  at  around  the  same  height.  This  suggests  that  it  is 
possible  that  the  ASSP  counts  are  not  liquid  droplets  but 
false  counts  from  light  scattered  off  the  snow  particles. 

In  order  to  answer  the  question  whether  the  ASSP 
detected  particles  are  supercooled  water  droplets.  the 
possible  sources  of  these  supercooled  droplets  need  be 
investigated.  There  are  two  possible  sources:  (1)  The  water 
droplets  being  initiated  below  the  freezing  level  and  then 
being  carried  upward  with  the  updraft  and  (2)  the  droplets 
being  initiated  in  situ  bg  supersaturation. 

In  order  to  investigate  the  possible  sources.  a  simple 
model  is  set  up.  Consider  a  closed  parcel  of  saturated  air 
initially  at  the  freezing  level  with  a  monodisperse 
distribution  of  liquid  water  content.  X-  •  and  the  initial 
radius  of  droplets  being  r  with  a  nucleus  of  NaCl  of  mass 
10-'*  gm.  In  this  parcel  of  air.  the  initial  snow  particle 
distribution  is  of  the  farm  No(h)e  ^,CdD.  The  upward 
velocity  of  the  parcel  is  taken  to  be  the  difference  between 
the  updraft  velocity  and  the  droplets  terminal  fallspeed. 
The  air  parcel  is  then  lifted  ad iabatical ly.  It  is  further 
assumed  that  there  is  no  entrainment  or  mixing  with  the 
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environmental  air.  The  output  is  the  radius  of  the 
supercooled  water  droplets*  vapor  saturation  ratio  with 
respect  to  ice  and  vapor  saturation  ratio  with  respect  to 
water.  According  to  Beard  and  Pruppacher  (1971),  the 
ventilation  factor  for  droplets  of  radius  20  micron  is  1.04 
and  so  the  ventilation  effect  has  been  neglected  in  the 
calculations. 


The  change  in  saturation  vapor  ratio  with  respect  to  ice 
is  computed  by  a  modified  form  of  Twomey,  1959. 


_  Lyi  ^  Lu  +  P  T _  Jty  j  €* 

Tree"  f  ~ 

-  +  _PikI 

R.T1  ct  Lit'e.  f. 


(3.  65) 


k 

6 

C 


■  supersaturation  with  respect  to  ice 

*  saturation  vapor  pressure  over  plain  ice 

*  saturation  vapor  pressure  over  plain  water 

■  temperature  of  air  parcel 

■  latent  heat  of  sublimation 

■  latent  heat  of  evaporation 


molecular  weight  of  water 
molecular  weight  of  air 


3  density  of  air 


-  gas  constant  of  air 

■  specific  heat  of  dry  air  at  constant  pressure 
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dH  *  change  of  height  of  air  parcel 
dW  ■  change  of  liquid  water  content 
dl  =  change  of  ice  content 

g  *  gravitational  acceleration  of  the  earth 
p  *  atmospheric  pressure 

The  first  term  on  the  right  hand  side  is  the  production  of 
water  vapor  due  to  lifting  and  cooling  of  the  air  parcel. 
The  second  term  is  the  depletion  (or  production)  of  water 
vapor  due  to  the  condensation  (or  evaporation)  of  water 
droplets.  The  third  term  is  the  depletion  of  water  vapor  due 
to  deposition  from  vapor  to  snow  particles. 


The  growth  of  liquid  droplets  is  computed  by  (Mason. 

1971) 


dS 


S<iU)  -  JjXtSif- 
ri&L _ -XI _ LlkJL 


-Lgw(_ir - 1)  +  -ft-J-1- 


(3.  66) 


where  r(h)  *  radius  of  water  droplets  as  a  funtion  of  height 
=*  density  of  water 
k  m  diffusivity  of  heat 

D  *  diffusivity  of  water  vapor 

ll  ■  number  of  ions  per  salt  molecule 
ms  =■  mass  of  salt 

M}  =  molecular  weight  of  salt 


dr  becomes  negative  when  S  is  negative  which  means  there  is 
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evpaoration.  The  change  of  liquid  water  content  is  given  by 


4  i  It 

-fll  r’  (hx>  -  —I} 


Assuming  the  snow  particles  are  circular  disks*  and  the 
mass  of  snow  is  related  to  its  diameter  by  x  *  o(  .  then 
the  growth  of  snow  particles  is  given  by  (Mason.  1971) 


_  1 )  +  2LE±.I- 

kT  pei 


(3. 67) 


Assuming  a  negative  exponential  distribution  of  snow 
particles,  the  change  in  ice  content  is  given  by 


dl 


1  cn U-fQot 


(3.  68) 


where  Ni>(  is  the  initial  intercept  of  snow-size  distribution 
and  N,xis  found  from 


loK  N.x-  lo».«  N<v.  ♦  <N.,  -N*  +dD>)v  (3.  69> 

with  N*  an  arbitrary  value.  The  temperature  of  the  air 
parcel  is  computed  by 

dT  *  ( - ) .  <3.  70) 

cf  cr 

The  terminal  fallspeed  of  the  droplets  is  assumed  to  be 
(Rogers.  1979): 


u  -  8000*r  (3.  71 ) 

where  u  *  droplet  terminal  fallspeed  in  cm/sec 
r  *  droplet  radius  in  cm 


The  upward  velocity  of  the  parcel  is  taken  to  be  the 
difference  between  the  updraft  velocity  and  the  terminal 
droplet  fallspeed. 


This  set  of  equations  is  solved  in  height  steps  and  the 
results  are  plotted  in  Figs.  3.  30  and  3.31  with  the  initial 
conditions  also  listed  in  the  figure  captions.  The  liquid 
water  content  used  is  1  gm/m*  and  is  extremely  high  for  this 
type  of  cloud.  A  monodisperse  cloud  with  droplets  of  radius 
20  microns  is  also  unrealistic.  In  real  life>  the  cloud 
droplet  distribution  is  approx imately  negative  exponential, 
which  means  that  there  are  more  small  drops  and  less  large 
drops.  Thus,  this  simple  calculation  simply  indicates  what 
happens  in  an  extreme  situation. 


Also  plotted  on  the  graphs  are  the  'water  saturation' 
lines  which  depict  the  amount  of  vapor  supersaturation  with 
respect  to  ice  if  the  vapor  pressure  is  kept  at  saturation 
with  respect  to  water.  The  results  show  that  the  vapor 
pressure  becomes  subsaturated  with  respect  to  water  starting 
a  few  decades  of  meters  above  the  freezing  level.  This  means 
that  supercooled  liquid  water  could  not  be  initiated  in  situ. 
For  small  updraft.  such  as  20  cm/s.  liquid  water  doplets 
cannot  be  maintained  for  a  long  distance  and  will  be 
evaporated  soon.  For  larger  updraft,  such  as  30  cm/s.  liquid 
water  droplets  will  be  maintained  up  to  1000  meters  above 


(«*»)  mBi»t 


0  10  20  30(Mm> 

Fig. 3. 30:  Change  of  saturation  vapor  pressure  iwth  respect  to  ice,  S  .  ,  satur¬ 
ation  vapor  pressure  with  respect  to  water,  S  and  superc6oled 
droplet  radius  with  height,  updraft  50  cm/s  ' 

Snow  particle  distribution  *  ^.le  10D 
Liquid  water  content  =  1  gm/m 
Initial  droplet  radius  =  20  microns 
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freezing  level.  Then  the  question  is  whether  the  droplet  can 
remain  being  supercooled  without  freezing  at  such 
temperature.  This  model  of  evporation  and 
condensation/deposition  permits  no  moisture  and  temperature 
mixing  between  the  parcel  and  the  environment  which  is  a 
reasonable  assumption  for  widespread  stratiform  cloud.  The 
results  from  this  simple  model  is  in  general  agreement  with 
complicated  models  such  as  Mason  and  Chien.  1962  and  Lee  and 
Pruppacher.  1977. 

In  any  case,  the  observed  occurrence  of  the  scant  liquid 
water  (Figs.  2.  8  -  2.11)  does  not  seem  to  correlate  with  the 
onset  of  the  collisional  stage.  Figs.  2.  12  and  2.13  clearly 
indicate  that  there  is  no  supercooled  liquid  water  for  the 
case  on  8  March.  1980.  However  there  is  still  the  third 
stage  of  snow  growth. 

Recent  studies  of  ice  particle  breakup  (Hobbs  and 
Farber.  1972 i  Vardiman.  1978)  emphasize  that  the  coexistence 
of  dense.  rapidly  falling  graupel  along  with  fragile 
aggregates  and  crystals  may  be  necessary  for  collisional 
breakup.  However,  the  observations  here  suggest  that  the 
presence  of  graupel  is  not  required  for  breakup.  The 
temperature  where  the  spectra  cease  evolution  is  about  -8*C 
which  is  in  the  temperature  range  of  scrolls  (Nakaya.  1954) 
and  hollow  prisms  (Mason.  1971).  It  could  be  that  these 


types  of  crystals  art  mors  fragile.  Martner  (1982)  also 
reported  the  occurrence  of  secondary  ice  crystal  production 
in  the  absence  of  supercooled  liquid  water. 


3.6  Oepos it ion— aggregation-rime  splintering  simulation 

In  order  to  study  the  possibility  of  the 
rime-splintering  mechanism  in  this  form  of  precipitation  an 
attempt  is  made  to  simulate  snow  growth  with  the  physical 
processes  of  deposition*  aggregation  and  rime-splintering. 


3.  6.  1  Two-parameter  formulation 

The  theoretical  representation  of  deposition  and 
aggregation  are  as  described  above.  According  to  the 
Hal iet-Mossop  theory*  the  rate  of  rime-splintering  depends  on 
the  environmental  liquid  water  content*  which  is  not 
available  in  the  theoretical  model.  Hence  rime-splintering 
is  represented  by  increasing  the  total  particle  concentration 
at  each  height  step  without  changing  the  total  mass.  The 
computational  procedure  is  that  first  the  changes  of  NQ  and 
are  computed  with  a  model  having  only  deposition  and 
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aggregation.  Where  IMtf  and  'h  arrives  at  their  observational 
'equilibrium'  positions!  the  total  particle  concentration  is 
made  to  increase  by  a  fixed  amount#  i.  e.  . 

- =  K  (3.  72) 

Ax  Aj 

In  order  to  keep  the  mass  constant#  the  change  in  N(  and  }\  is 
given  by: 

N,  -  N.  (3.73) 

Xx  =■  C<~  +  ~->  X?*  (3.74) 

Ni  Ai 

In  the  computations,  the  values  of  K  used  are  based  on 
the  splintering  production  rates  listed  in  Table  3.  3 
multiplied  by  the  time  elapsed  in  one  height  step.  The 
results  for  these  computations  for  the  three  ASD  are  plotted 
in  Figs.  3.  32#  3.33  and  3.34.  The  value  of  physical 

parameters  are  listed  in  Table  3.6.  It  shows  the  increase  in 
total  concentration  for  certain  initial  N,  and  \  and  change  of 
concentration.  The  results  from  the  model  indicate  that 
riming  can  stop  the  decrease  of  X  for  a  certain  while  with 
Na increasing.  Then  both  N.  and  \  decrease  unceasingly.  This 
formulation  fails  to  predict  the  equilibrium  X  that  is 
observed.  This  is  a  further  indication  that  the  third  stage 
of  snow  growth  in  the  observat ional  data  cannot  be  caused  by 
rime-splintering. 


>arameter  model  of  depi 
ig,  26  February,  1980 
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Table  3.6:  Values  of  parameters  used  in  aggregation-deposition- 
rime  splintering  simulations  (c.g.s.  units) 
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3.  6.  2  Three-parameter  formulation 


In  a  three-parameter  formulation*  three  moment 
conservation  equations  are  derived.  This  means  that  the 
changes  of  the  total  particle  concentration*  the  total  mass 
and  the  total  radar  reflectivity  factor  are  being  completely 
controlled  by  the  model.  This  makes  the  three-parameter 
formulation  ideal  to  study  the  effect  of  rime-splintering* 
which  affects  only  the  total  particle  concentration  without 
affecting  the  total  mass  and  the  total  radar  reflectivity 
factor.  A  three-parameter  formulation  with  the  physical 
processes  of  vapor  deposition*  aggregation  and 
rime-splintering  results  in  the  following  equations. 


The  change  of  particle  concentration  flux  is  given  by 


-  -  iiLilA./  f  x'-y*'  <  x+v  j1  ;  xb -yl 

1  i 


I  .-"V  dXdY 


E.  T<i no  r ((*«») 

•\  k+Vt3 


(3.  75 > 


where  the  first  term  on  the  right  hand  side  is  the  depletion 
of  particles  due  to  aggregation  and  the  second  term  is  from 
(3.61)  and  is  the  production  of  particles  due  to 
rime-splintering.  The  other  two  moment  conservation 
equations  are  similar  to  those  presented  earlier  except  that 
collisional  breakup  is  being  ignored.  The  derivation  results 
in  two  implicit  equations  in  ^»and  <T  * 


(I  rc\>*  <r+0  -i 
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where  I,  end  Ij  ere  from  (3.29)  end  (3.31). 

In  the  computations*  the  velues  of  n'.  the  concentretion 
of  droplets*  is  teken  from  the  ASSP  concentretion  of 
perticles  greeter  then  24  microns.  The  velues  of  the  other 
peremeters  ere  identicel  to  those  listed  in  Teble  3.4.  The 
results  using  the  peremeter  velues  of  the  29  Feb.  *  1980  case 
ere  plotted  in  Figs.  3. 39-3. 37.  The  computed  evolution  of  N^. 
X*end  <r  do  not  metch  the  observetions  desirably*  especially 
for  O'  .  The  predicted  <r  is  first  negative*  then  positive  and 
then  negative  again*  while  the  observational  value  of  <T  is 
first  negative  end  then  positive. 

In  the  computations  of  this  deposition,  aggregation  and 
rime-splintering  model  using  the  condition  of  the  8  March. 
1980  case*  it  is  found  that  \  will  become  close  to  zero. 
When  it  does  that*  the  computations  become  numerically 
unstable.  In  order  to  circumvent  this  problem  the 
calculations  are  performed  using  a  deposition*  aggregation* 
collisional  breakup  and  rime-splintering.  The  model  is 
computed  twice.  The  first  time  uses  the  observed  value  of 
n'*  the  number  density  of  droplets  greater  than  24  microns 
and  a  small  value  of  the  breakup  parameter*  c  -  0.009 
(Figs.  3.  38*  3.40  and  3.  42).  The  reason  for  using  a  small 
value  of  c  is  that  on  the  one  hand*  the  effect  of  breakup 
will  not  overshadow  the  effect  of  rime-splintering*  while  on 
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Fig. 3. 38:  N*  from  three-parameter  model  of  deposition,  aggregation  and  rime 
splintering  with  weak  breakup,  8  March,  1980 


Fig. 3. 39:  Similar  to  Fig. 3. 38  but  with  no  rime-splintering 
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the  other  hand. 

the 

effect 

of  breakup 

will 

prevent  A* 

from 

becoming  too  small. 

thus 

causing  numerical 

instability. 

The 

second  computation 

is 

done 

with  n'  = 

0 

and  c  =  0 

.  005 

(Figs.  3.  39.  3.41 

and 

3.  43). 

The  results 

from  the 

two 

computation*  are  then  compared  with  each  other  and  with  the 
observations.  By  doing  this  comparison!  the  effect  of 
rime-splintering  can  be  seen.  The  results  from  the  two 
computations  show  that  rime-splintering  has  a  slight  effect  on 
the  evolution  of  At  and  <T  but  the  effect  is  not  adequate  to 
explain  the  evolution  depicted  in  the  observations. 


3.  7  Effects  of  melting 

The  evolution  of  snow-size  spectra  to  rain-size  spectra 
due  to  melting  can  be  studied  by  assuming  that  all  moment 
fluxes  are  conserved  during  melting.  This  is  equivalent  to 
saying  that  the  term*  on  the  right  hand  side  of  (3.1)  are 
identically  zero.  For  a  two-parameter  formulation.  the 
conservation  of  mass  flux  is  represented  by. 

JatDk,o(D?  Nse*AtI>  dD  -  JmT  y  f  D1  Nre'A?I>  dD  (3.78) 

where  the  subscript  s  refers  to  snow  and  the  subsript  r 
refers  to  rain  and  a.  b.  and  ^  are  the  same  notation*  as  in 
(3.20)  and  (3.21).  The  conservation  of  radar  reflectivity 
flux  is  represented  by. 
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Among  the  data  from  the  spirals*  there  is  only  one  case 
in  which  observations  were  made  where  snow  had  been 
completely  melted.  The  data  in  this  case*  which  is  the 
second  spiral  on  17  April*  1981*  is  used  to  test  the  effect 
of  melting  on  the  spectrum.  Using  (Locatelli  and  Hobbs* 
1974*  Rogers*  1979  and  in  c.  g.  s.  units)  a^  =  100* 
bs  -  O.  IS.  ar  -  1421*  h,.  *  0.  5.  *  -  0.  002.  JS  »  2,  the 
theoretically  derived  values  for  the  melted  spectrum  can  be 
computed.  Fig.  3. 44  shows  the  two-parameter  fit  to  the 
observed  spectra  before  and  after  melting.  The  spectrum 
predicted  from  theory  is  also  plotted  for  comparison.  The 
slope  of  the  observed  melted  spectrum  is  shallower  than  the 
slope  from  theory  while  the  intercepts  are  almost  the  same. 
The  observed  mass  flux  is  about  twice  that  of  the  computed 
mass  flux.  By  varying  the  values  of  the  parameters.  it  is 
found  that  the  mass  flux  is  rather  sensitive  to  at»  the 
fallspeed  coefficient.  If  a^  ■  290  is  used*  the  computed 
mass  flux  is  more  comparable  to  the  observed  value.  This 
might  suggest  that  the  terminal  fallspeed  of  the  type  of  snow 
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immediately  above  the  freezing  level  is  larger.  The  computed 
spectrum  from  using  a$  *  250  is  also  plotted  on  Fig.  3. 44. 
Comparing  this  spectrum  uiith  the  observed  one>  it  can  be 
deduced  that  during  melting,  aggregation  is  more  efficient 
than  breakup.  Thus  small  drops  are  depleted  and  large  drops 
formed. 


For  a  three-parameter  formulation>  the  conservations  of 
concentration  flux>  mass  flux  and  radar  reflectivity  flux  are 
given  by. 
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Then  Ny*  \  and  $"T  are  obtained  by, 
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Fig.  3.  43  show*  the  three-parameter  spectra  before  and 
after  melting.  The  results  are  similar  to  those  obtained 
from  the  two-parameter  formulations.  Again  the  physical 
implications  from  these  results  are  similar  to  those  obtained 
from  the  two-parameter  analyses.  It  is  interesting  to  see 
that  during  melting<  the  resulting  raindrop  size  spectrum  is 
always  subexponential. 


-  parameter  melting 


Fig. 3. 45:  Three-parameter  fit  to  melting 
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CHAPTER  4 

CONCLUSIONS 


Both  the  observational  data  and  the  theoretical  models 
indicate  clearly  that  snow  growth  goes  through  three  distinct 
processes*  vapor  deposition  stage*  aggregation  stage  and  a 
secondary  production  stage.  Initially  when  most  of  the  snow 
particles  are  small*  they  tend  to  grow  by  water  vapor 
depositing  directly  onto  them.  By  the  time  some  of  the 
particles  grow  to  large  enough  sizes*  thus  increasing  the 
collisional  probability*  aggregation  sets  in.  Aggregation 
causes  the  depletion  of  small  particles  and  the  increase  in 
large  particles.  The  average  size  of  snow  particles 
therefore  increases  rapidly. 

The  third  stage  of  snow  growth  is  the  secondary 
production  stage.  Thera  has  been  much  discussion  on  what 
causes  the  secondary  production.  There  are  two  hypotheses* 
namely  rime-splintering  and  collisional  breakup.  An 
important  ingredient  of  rime-splintering  is  supercooled 
liquid  water.  The  observational  data  indicate  that 
supercooled  liquied  water  is  present  in  two  of  the  three 


flights  considered  and  is  absent  in  the  other  flight.  But* 
the  secondary  production  stage  is  still  observed  in  all  three 
flights.  Moreover<  in  the  flights  when  supercooled  liquid 
water  is  observed!  the  occurrence  of  supercooled  liquid  water 
is  rare  and  apparently  not  correlated  with  any  particular 
stage  of  snow  growth.  Since  the  secondary  production  stage 
can  occur  even  in  the  absence  of  supercooled  liquid  water! 
rime-splintering  is  unlikely  to  be  a  significant  candidate 
for  causing  secondary  production.  On  the  other  hand!  the 
fact  that  the  onset  of  the  secondary  production  stage  causes 
the  snow-size  spectra  to  have  an  equilibrium  slope  is 
indicative  of  collisional  breakup.  Physically!  aggregation 
of  snow  particles  increases  the  sizes  of  the  particles.  Buti 
when  the  particles  become  too  big(  they  become  more  fragile. 
Thus  they  tend  to  break  up  upon  collisions. 


Theoretical  investigations  produce  similar  conclusions. 
Theoretical  models  with  the  processes  of  vapor  deposition! 
aggregation  and  collisional  breakup  give  results  comparable 
to  the  observational  data!  while  theoretical  models  with  the 
processes  of  vapor  deposition!  aggregation  and 
rime-splintering  fail  to  give  results  similar  to  the  data. 
Therefore!  both  theory  and  observation  point  out  the 
importance  of  collisional  breakup  in  snow  growth.  From  the 
formulations  of  the  collisional  breakup  term  in  the  theory 
and  the  results  of  the  laboratory  experiment!  it  can  be 
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Concluded  that  when  two  snow  particles  collide*  they  are  more 
likely  to  aggregate  than  to  break  up.  However*  when  they  do 
break  up*  they  generate  fragments  distributed  in  a  somewhat 
negative  exponential  manner*  the  larger  the  parent 
particles*  the  more  numerous  the  fragments.  Another 
important  finding  in  the  theoretical  studies  is  that  crystal 
types  need  not  be  taken  into  account  in  order  to  produce  the 
different  stages  of  snow  growth*  although  from  the  2-D  data 
it  is  evident  that  different  crystal  types  dominate  at 
different  stages. 

Data  from  two  of  the  spirals  indicate  that  when  the 
atmosphere  is  subsaturated*  evaporation  of  the  particles  will 
prohibit  the  occurrence  of  an  equilibrium  distribution  slope. 
A  possible  explanation  is  that  the  number  of  large  particles 
is  being  maintained  by  aggregation*  while  collisional  breakup 
generates  minute  particles  that  are  not  detectable  by  the  1-D 
precipitation  probe.  Under  conditions  of  deposition*  these 
small  particles  will  grow  into  detectable  sizes*  but*  under 
conditions  of  subsaturation*  these  small  particles  not  only 
cannot  grow  into  detectable  sizes*  but  rather  they  are 
evaporated.  This  causes  the  spectrum  to  become  'flatter'  and 
so  the  slope  will  decrease  below  the  10  cm1  limit  observed  in 
the  other  spirals.  Since  only  the  data  from  two  spirals* 
both  taken  on  the  same  day*  are  available,  this  finding  can 
be  considered  to  be  preliminary.  More  investigations  of  snow 
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growth  in  subsaturated  conditions  arc  needed. 

The  theoretical  models  show  that  increasing  the  initial 
concentration  without  increasing  the  initial  mass  (physically 
equivalent  to  seeding)  causes  the  onset  of  the  aggregation 
stage  of  snow  growth  to  start  slightly  earlier  but  delays  the 
onset  of  the  third  stage  by  a  more  significant  amount. 
Physically  the  increase  in  concentration  without  increasing 
the  mass  means  that  the  number  of  collisions  increases. 
Therefore/  the  aggregation  stage  starts  earlier.  On  the 
other  hand/  increasing  the  concentration  without  increasing 
the  mass  means  that  the  average  size  of  particles  is  smaller. 
So/  it  would  take  longer  for  them  to  grow  to  such  a  size  that 
collisional  breakup  becomes  efficient.  Since  neither 
aggregation  nor  collisional  breakup  affects  the  total  mass/ 
so  the  final  slope  does  not  change. 

The  three-parameter  negative  exponential/  power  law  size 
distribution  provides  a  more  sophisticated  but  yet 
analytically  feasible  method  of  describing  snow  size  spectra. 
This  new  formulation  enables  the  investigator  to  have  one 
more  degree  of  freedom  in  the  analysis.  The  difficulty  in 
using  this  formulation  lies  in  the  instrumentation  because 
there  is  a  great  discrepancy  between  the  particle 
concentration  measured  by  the  1-D  cloud  probe  and  the  1-D 
precipitation  probe.  A  three-parameter  analysis  of  data  is 
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very  sensitive  to  the  smell  end  of  particle  distribution. 
When  the  discrepancy  between  the  cloud  probe  and  the 
precipitation  probe  is  resolved.  the  three-parameter 
formulation  will  become  very  useful. 

When  shown  in  the  N0-)\  space.  it  is  seen  that  at  the 
third  stage  of  snow  growth.  the  different  cases  all  have  the 
same  slope  but  different  intercepts.  The  differences  in 
intercept  are  obviously  tied  to  the  precipitation  rate  and 
are  possibly  due  to  the  different  dynamical  conditions. 

The  studies  here  make  no  attempt  to  explain  the  initial 
nucleation  of  snow  particles.  For  instance,  the  theoretical 
model  assumes  a  certain  snow-size  distribution  as  initial 
condition.  A  full  model  of  snow  growth  should  really  start 
out  with  primary  nucleation.  Unfortunately.  the  current 
state  of  ice  ncleation  is  still  not  adequate  to  provide 
enough  insight  to  incorporate  this  into  these  models. 

It  would  be  good  if  different  investigators  could  obtain 
similar  data  from  various  types  of  winter-storms  so  as  to 
provide  a  climatology  of  snow-gowth  processes.  This  could 
facilitate  the  study  of  the  microphysical  processes 
associated  with  different  storm  dynamics.  Future  studies  of 
snow  growth  processes  can  be  pursued  by  using  the  2-D  data 
quantitatively.  Right  now.  investigators  who  use  the  2-D 


data  tand  to  reduce  the  2-D  data  to  1-D.  Perhaps*  another 
approach  is  to  express  snow-size  distribution  as  a  function 
of  area*  f(A)dA<  instead  of  the  conventional  function  of 
diameter*  f(D)dD. 
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ABSTRACT 

In  studies  of  precipitation  growth,  comparisons  between  theory  and  observation  are  difficult  because  of 
the  problem  in  obtaining  a  complete  4-dimensional  (space  and  time)  description  of  the  kinematic,  ther¬ 
modynamic  and  microphysical  properties  of  the  atmosphere.  A  new  flight  plan  has  been  devised  which 
permits  one  to  observe  the  height  evolution  of  snow-size  spectra  in  a  reference  frame  where  the  effects  of 
horizontal  gradients  and  temporal  changes  are  minimized.  The  flight  plan,  termed  the  advecting  spiral 
descent  (ASD),  requires  an  aircraft  to  start  aloft  in  a  mesoscale  precipitation  area  and  then  spiral  downward 
in  a  constant  bank  angle,  descending  at  approximately  the  mean  fallspeed  of  snow.  Qualitative  comparisons 
between  ASD  observations  and  particle  growth  theory  suggest  that  snow  evolves  through  at  least  three 
stages  characterized  by  deposition,  aggregation  and  breakup.  The  breakup  process  serves  to  limit  the  number 
of  large  snow  particles  and  interacts  with  aggregation  to  produce  a  limiting  value  of  the  slope  of  the  snow- 
size  spectrum. 


1.  Introduction 

The  ultimate  test  of  a  theoretical  model  of  a  phys¬ 
ical  process  is  how  well  the  model  can  describe  the 
natural  occurrence  of  the  process.  In  studies  of  pre¬ 
cipitation  growth,  such  comparisons  are  difficult  be¬ 
cause  of  the  problem  in  obtaining  a  complete  4-di¬ 
mensional  (space  and  time)  description  of  the 
kinematic,  thermodynamic  and  microphysical  prop¬ 
erties  of  the  atmosphere.  The  rather  detailed  de¬ 
scriptions  of  the  atmosphere  obtained  from  micro¬ 
physical/dynamical  models  are  typically  beyond  our 
observational  capability  for  verification. 

The  theoretical  modeling  of  snow  growth  by  one 
of  the  authors  (Passarelli,  1978a,b)  raised  the  ques¬ 
tion  of  how  to  verify  the  model  results  using  mea¬ 
surements  of  snow-size  spectra  obtained  by  a  single 
aircraft  equipped  with  laser  imaging  probes.  In  this 
paper  we  present  a  flight  plan  which  permits  one  to 
observe  the  height  evolution  of  snow-size  spectra  in 
a  reference  frame  where  the  effects  of  horizontal 
gradients  and  temporal  changes  are  minimized.  The 
aircraft  is  constrained  to  drift  with  a  region  of  falling 
snow  in  a  Lagrangian  frame  of  reference  tied  to  the 
average  particle  motion.  The  height  evolution  of  the 
size  distribution  of  snow  then  reduces  to  a  simple  1- 
dimensional  steady-state  problem  so  that  a  single 
aircraft  can  provide  sufficient  data  for  compari¬ 
son  with  results  from  1 -dimensional  microphysical 
models. 

First,  we  present  a  simple  sedimentation  model  of 
falling  snow  to  provide  the  physical  rationale  behind 
the  flight  plan,  followed  by  a  description  of  the  air¬ 
craft  flight  track.  Several  examples  of  data  are  given 


along  with  a  physical  interpretation  of  the  results. 
A  detailed  quantitative  comparison  between  theory 
and  observation  will  be  presented  in  a  later  paper. 

The  results  show  that  snow  growth  is  a  well-be¬ 
haved  phenomenon  whose  evolutionary  processes 
can,  to  a  large  extent,  be  observed  with  current  in¬ 
strumentation.  Snow  growth  is  apparently  dominated 
initially  by  deposition,  followed  by  aggregation  and 
then  breakup.  The  results  suggest  that  an  aggrega¬ 
tion-breakup  equilibrium  is  acting  to  produce  a  lim¬ 
iting  value  of  the  slope  of  the  size  distribution  of 
snow. 

2.  The  physical  rationale  behind  the  sampling  pro¬ 
cedure 

The  change  in  the  size  distribution  of  snow  at  any 
point  in  time  and  space  is  governed  by  advection  (air 
motion  and  particle  fallspeed),  primary  particle 
growth  (deposition  or  evaporation  and  riming  of  su¬ 
percooled  drops),  the  redistribution  of  mass  within 
the  size  distribution  (aggregation  and  breakup)  and 
the  production  of  new  ice  particles  (nucleation).  The 
question  addressed  here  is  whether  a  single  aircraft 
can  sample  the  atmosphere  in  such  a  way  so  as  to 
separate  the  effects  of  advection  from  these  micro- 
physical  processes. 

Perhaps  the  most  desirable  situation  for  a  single 
aircraft  study  of  snow  growth  is  that  of  widespread, 
steady-state  snow.  In  the  absence  of  horizontal  gra¬ 
dients,  horizontal  advective  processes  make  no  con¬ 
tribution  to  changes  in  the  size  distribution.  An  air¬ 
craft  can  merely  sample  the  atmosphere  at  various 
heights,  at  leisure,  provided  that  conditions  are 


0022-4928/82/040697- 10S06. 50 
©  1982  American  Meteorological  Society 


The  u.S.  Government  1s  authorized  to  reproduce  and  tell  this  report. 
Permission  for  further  reproduction  by  others  must  be  obtained  from 
the  copyright  owner. 


698 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  39 


z 


Fig.  I.  Particle  trajectories  from  steady-state  line  source,  scales  in  km. 


steady.  The  particle  growth  can  then  be  inferred  from  the  large  ( 1 30  cm  s-')  particles  is  indicated  by  hor- 
the  height  change  in  the  size  distribution.  Unfortu-  izontal  hatching  and  the  region  occupied  by  the  small 
nately,  even  in  widespread  storms  there  are  usually  (80  cm  s"')  particles  is  indicated  by  vertical  hatching, 
mesoscale  areas  of  precipitation  which  move  and  The  right-hand  side  of  this  precipitation  trail  con- 
change  with  time.  Horizontal  gradients  in  the  prop-  tains  no  small  particles  while  the  left-hand  side  con- 

erties  of  the  atmosphere  coupled  with  the  effects  of  tains  no  large  particles.  However,  the  region  of  over- 

advection,  fallspeed  dispersion  and  wind  shear  create  lap  for  these  two  fallspeeds  has  uniform  properties 
a  difficult  sampling  problem  even  if  conditions  are  which  are  identical  to  the  source  region.  If  the  par- 
steady.  tide  fallspeeds  are  bounded  by  these  hypothetical 

To  illustrate  this,  consider  a  simple  two-dimen-  limits  then  the  complete  size  distribution  (all  fall- 
sional  sedimentation  model  of  a  finite,  horizontally  speeds)  will  be  constant  in  the  region  of  overlap.  The 
homogeneous,  steady,  line  source  of  snow  aloft  which  properties  of  the  snow  in  the  central  region  in  Fig. 
is  10  km  wide  and  located  5  km  above  the  earth’s  1  are  identical  to  those  corresponding  to  an  infinite 
surface  in  an  atmosphere  having  constant  shear  (e.g.,  line  source  aloft. 

Marshall,  1953).  We  assume  that  no  particle  growth  This  example  illustrates  that  for  a  steady,  finite, 
occurs  and  that  updrafts  are  weak  compared  to  upper-level  source  region  in  which  the  dynamical  and 
fallspeeds.  A  particle  starting  from  the  source  at  microphysical  properties  are  relatively  constant  in 
(Wight  (Z0)  will  fall  and  be  deflected  from  its  initial  horizontal  space,  there  will  be  a  limited  region  below 
h&rizontal  position  (X0)  with  respect  to  the  source,  that  behaves  as  if  there  were  an  infinite,  horizontally- 
i.e.,  the  particle  trajectory  is  homogeneous  source  aloft.  In  this  region  the  effects 

2  of  horizontal  advection  can  be  ignored  and  snow 

X*.  x  +  -  (Z°  ~  (1)  growth  can  be  treated  as  a  1-dimensional  (height), 

0  .2  Vf  '  1  steady-state  problem.  Of  course  the  smaller  the  hor¬ 

izontal  gradients,  the  greater  the  depth  over  which 
where  X  is  the  horizontal  position  of  the  particle  with  this  assumption  will  remain  valid, 
respect  to  the  source,  S  the  shear,  Z  the  height  below  The  problem  of  temporal  variability  of  the  source 
the  source  and  Vf\Yic  fallspeed.  Neglecting  horizon-  region  can  be  examined  in  a  similar  manner.  Let  us 
tal  gradients  in  the  wind,  all  particles  having  the  again  consider  only  two  fallspeeds,  80  and  130  cm 
same  fallspeed  will  have  parallel  trajectories.  Fig.  1  s"‘,  and  use  the  identical  geometry  and  shear  as  in 

is  a  scale  drawing  of  the  precipitation  “trails”  for  the  previous  case.  However,  let  the  steady  source  act 
two  particle  fallspeeds  (80  cm  s~'  and  130  cm  s~')  for  only  3125  s.  During  this  time,  the  small  flakes 
assuming  a  constant  wind  shear  of  10*J  s“\  Note  fall  2.5  km  and  the  large  flakes  fall  4.0  km.  Fig.  2 
that  the  reference  frame  is  fixed  to  the  source.  The  shows  the  regions  occupied  by  the  two  populations 
width  of  the  region  occupied  by  particles  of  a  given  after  the  source  terminates  at  3125  s,  and  at  6250 
fallspeed  is  always  equal  to  the  width  at  the  source,  s  when  the  first  small  flake  hits  the  ground.  After 
However,  because  particles  of  different  fallspeeds  3125s  there  is  a  large  region  in  which  the  snow  has 
define  different  trajectories,  size  sorting  occurs  and  the  steady-state  properties  of  the  initial  source.  At 
horizontal  gradients  develop.  The  region  occupied  by  6250  s,  31 25  s  after  the  source  is  stopped,  there  is 


Fig.  2.  Particle  trajectories  from  time-dependent  line  source.  Kales  in  km. 
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still  a  region  close  to  the  surface  which  has  the 
steady-state  properties  of  the  source.  Recall  that  this 
region  also  has  the  properties  of  a  horizontally  in¬ 
finite  source. 

This  example  illustrates  that  for  a  source  which 
varies  slowly  in  time,  there  can  be  a  region  below 
that  behaves  as  if  the  source  were  infinite  and  steady. 
The  center  of  this  region  falls  at  a  rate  equal  to  the 
average  particle  fallspeed.  Ideally  if  one  could  use 
an  aircraft  to  follow  this  region,  then  the  snow  growth 
problem  reduces  to  a  1 -dimensional  (height)  steady- 
state  problem.  We  will  refer  to  this  region  as  the  1- 
dimensional  region. 

The  primary  assumption  in  both  of  these  analyses 
is  that  the  fall-speeds  lie  between  finite  limits.  Also 
the  lower  limit  cannot  be  zero.  These  are  fairly  good 
assumptions  for  snow  since  observations  show  that 
the  fallspeed  of  snow  is  only  a  weak  function  of  par¬ 
ticle  size.  This  is  particularly  true  of  unrimed  crystals 
and  aggregates  (e.g.,  Locatelli  and  Hobbs,  1974). 

To  illustrate  the  narrow  dispersion  of  snowflake 
fallspeeds  we  can  consider  a  specific  example  of  an 
exponential  distribution  of  equivalent  melted  diam¬ 
eters  and  a  power  law  f allspeed-melted  diameter  re¬ 
lationship,  i.e., 

N(Dm)  =  Ntf-™-,  (2) 

Vf=aDm\  (3) 

where  N(D„)M)„  is  the  concentration  of  particles 
in  the  melted  diameter  interval  [Dm,  Dm  +  A Dm\,  N0 
and  X  are  the  distribution  parameters,  Vf  is  the 
fallspeed  and  a  and  b  depend  on  the  snowfall  type. 
One  can  show  that  the  relative  dispersion  (standard 
deviation/mean)  of  the  mass  flux  as  a  function  of 
fallspeed  depends  on  b  only,  i.e., 

<r  f~r(4  +  3fi)r(4  4-  b)  l,/a 

V  L  r*(4  +  2A)  J  ’  v  ; 

where  a  is  the  standard  deviation  of  the  mass  flux 
versus  fallspeed  distribution,  V  is  the  mean  flux 
fallspeed  and  r  represents  the  gamma  function.  For 
aggregate  snow  A  0.3  (Langleben,  1934),  so  that 
the  relative  dispersion  is  only  1 5%.  In  the  sedimen¬ 
tation  examples,  we  used  80  and  130  cm  s~'  as  the 
fallspeed  limits.  The  ±2<r  fallspeeds  are  74  and  1 37 
cm  s-'  assuming  a  mean  mass  flux  fallspeed  of  103 
cm  s-1  and  b  =»  0.3.  Thus  for  snow,  the  vast  majority 
of  the  mass  flux  is  contained  within  a  narrow  range 
of  fallspeeds. 

Summarizing  the  discussion  of  the  simple  sedi¬ 
mentation  models  we  can  say  that  there  is  a  region 
below  a  time-dependent,  horizontally  finite  source 
which  behaves  as  if  the  source  were  steady  and  in¬ 
finite.  The  sampling  scheme  suggested  by  these  anal¬ 
yses  is  to  start  at  the  center  of  such  a  source  region 
and  then  follow  the  motion  of  an  average"  particle 
as  it  falls  and  is  advected  by  the  horizontal  wind. 


Observations  of  the  particle-size  distribution  in  this 
reference  frame  can  then  be  interpreted  by  means 
of  a  1-dimensional  (height),  steady-state  model  of 
particle  growth. 

The  technique  is  likely  to  succeed  in  regions  of 
large-scale  winter  storms  in  which  vertical  air  veloc¬ 
ities,  horizontal  gradients  and  temporal  variations 
are  small.  If  these  conditions  prevail  then  there  is 
usually  very  little,  if  any  supercooled  water  so  that 
riming  growth  is  not  significant  (e.g.,  Passarelli, 
1978b;  Herzegh  and  Hobbs,  1980).  The  precipitation 
is  in  the  form  of  snow  and  ice  crystals,  and  nucle- 
ation,  deposition,  aggregation  and  possibly  breakup 
are  the  microphysical  processes  which  can  act  to 
change  the  size  distribution. 

3.  The  aircraft  flight  track 
a.  Flight  procedure 

Having  defined  an  optimal  region  for  studying 
snow  growth  in  the  natural  atmosphere,  we  propose 
the  following  flight  profile  for  sampling  with  a  single 
aircraft.  Starting  aloft  in  a  mesoscale  precipitation 
area,  the  aircraft  is  placed  in  a  constant  bank  angle 
(~15°)  and  a  constant  descent  rate  (—  1  m  s_l). 
The  aircraft  spirals  dowr  yard  at  approximately  the 
mean  fallspeed  of  snow  and  the  loops  of  the  spiral 
drift  with  the  wind.  Ideally,  if  conditions  are  quasi¬ 
steady  and  the  properties  of  the  atmosphere  are  fairly 
uniform  over  a  length  scale  somewhat  larger  than 
the  diameter  of  the  loops,  then  the  aircraft  will  re¬ 
main  in  a  l -dimensional  region.  We  call  this  flight 
track  the  advecting  spiral  descent  (ASD). 

An  important  advantage  of  the  ASD  in  sampling 
a  I -dimensional  region  is  that  one  does  not  require 
an  a  priori  knowledge  of  the  vertical  wind  profile. 
The  aircraft  adjusts  automatically  to  the  wind  shear 
since  the  loops  of  the  spiral  are  flown  relative  to  the 
air  (bank  angle  constant)  rather  than  relative  to  the 
ground.  In  principle,  knowing  the  wind  as  a  function 
of  height  and  an  average  fallspeed  one  could  con¬ 
struct  a  series  of  horizontal  passes  which  would  ac¬ 
complish  the  same  task,  but  the  procedure  is  rela¬ 
tively  complex  in  four-dimensions  as  compared  to  the 
ASD  approach. 

The  analysis  of  particle  size  spectra  sampled  via 
the  ASD  technique  can  be  performed  by  averaging 
spectra  over  a  complete  loop  of  the  spiral.  A  loop  is 
defined  as  when  the  aircraft  completes  a  360°  turn. 
This  serves  to  average  any  horizontal  inhomogene¬ 
ities.  Ideally,  one  would  like  horizontal  inhomogene¬ 
ities  to  be  minimal  over  the  radius  of  the  loops. 

Another  approach  is  to  compare  particle  size  spec¬ 
tra  at  various  heights  which  occurred  in  the  same 
sector  of  different  loops  (e.g.,  using  magnetic  head¬ 
ing).  Ideally,  if  all  the  loops  are  the  same  size  and 
the  aircraft  performs  a  1  ms'1  ASD,  then  each  point 
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Angle  (Degrees) 

Fio.  3.  Turn  radius  as  a  function  of  bank  angle  for  various  true  air  speeds.  Air  speeds  in  m  s'1. 


of  the  aircraft  trajectory  corresponds  to  the  trajec¬ 
tory  of  a  1  m  s~*  particle.  After  turning  360°  (com¬ 
pleting  the  next  loop)  the  aircraft  should  encounter 
the  same  1  m  s~‘  particle  at  a  lower  altitude.  This 
is  stretching  the  technique  to  the  ultimate  limit.  In 
practice,  as  discussed  in  the  next  section,  the  loops 
get  smaller  as  the  aircraft  descends  so  that  the  air¬ 
craft  does  not  actually  re-encounter  the  same  1  m 
s'1  particle.  Nevertheless,  the  observed  point-to-point 
correlations  can  be  quite  high  as  we  shall  show. 


b.  The  ASD  flight  characteristics 


During  a  coordinated  turn  at  a  constant  bank  an¬ 
gle,  the  aerodynamic  lift,  gravitational  and  centrif¬ 
ugal  forces  are  balanced.  One  can  show  that  the  ra¬ 
dius  of  a  turn,  r  is 


r 


(TAS)2 

g 


cold. 


(5) 


where  TAS  is  the  aircraft  true  air  speed,  9  is  the 
bank  angle  and  g  is  the  gravitational  acceleration. 
The  turn  radius  as  a  function  of  bank  angle  for  var¬ 
ious  true  air  speeds  is  plotted  in  Fig.  3.  During  a 
descent,  a  pilot  will  usually  maintain  a  constant  in¬ 
dicated  air  speed  (IAS)  which  is  related  to  the  true 
air  speed  by 

Po(IAS)J  =  p(TAS)2,  (6) 


where  p0  is  the  density  of  air  at  a  standard  temper¬ 


ature  and  pressure  (0°C  and  1013  mb)  and  p  is  the 
actual  air  density.  In  a  descent  where  the  bank  angle 
and  indicated  air  speed  remain  constant,  the  true  air 
speed  is  less  at  lower  levels  and  the  radius  of  the 
loops  decreases  in  accordance  with 


r 


Po  (IAS)2 
P  g 


cottf. 


(7) 


Hence,  if  p  doubles  over  the  depth  of  the  spiral,  the 
turn  radius  will  decrease  by  one-half. 

Ideally,  one  would  like  to  have  the  turn  radius  as 
small  as  possible  to  minimize  the  effects  of  horizontal 
in  homogeneities.  Practically,  the  turn  radius  is  dic¬ 
tated  by  the  pilot’s  ability  to  safely  maintain  a  given 
bank  angle  for  long  periods  of  time  in  instrument 
flight  conditions.  A  bank  angle  of  30s  is  probably  a 
maximum  upper  limit.  In  the  examples  which  follow, 
we  used  bank  angles  of  15s  to  20°. 


4.  Examples  of  ASD  data 

a.  Data  and  Instrumentation 

To  illustrate  the  technique  we  present  data  from 
the  Air  Force  Geophysics  Laboratory  instrumented 
C-l  30  cloud  physics  aircraft  which  is  equipped  with 
a  variety  of  sensors.  Only  the  data  from  the  Particle 
Measuring  Systems  (PMS)  200-Y  probe  will  be  dis¬ 
cussed  here.  This  is  a  1 -dimensional  laser  imaging 
probe  which  counts  and  sizes  particles  into  IS  size 
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Table  I.  Summary  of  spirals. 


Storm 

Spiral  top 

Spiral  base 

Height 

Tempera¬ 
ture  CC) 

Height 

Tempera¬ 
ture  co 

1 

6740  m 

-29 

3060  m 

-6 

2 

6940  m 

-26 

2320  m 

0 

3 

5860  m 

-21 

3630  m 

-5 

categories  ranging  from  300  to  4500  nm,  each  300 
pm  wide.  Three  flights  are  presented.  Two  were  made 
off  the  coast  of  Washington  (25  and  26  February 
1980)  and  one  off  the  coast  of  New  Hampshire  (8 
March  1980).  One  ASD  was  performed  on  each  day 
and  these  will  be  referred  to  as  spirals  1,  2  and  3 
respectively. 

All  three  spirals  took  place  in  winter  cyclonic 
storms.  Spiral  1  was  flown  just  ahead  of  an  occluded 
front  and  spiral  2  was  flown  ~  100  km  east  of  a  low- 
pressure  center.  Spiral  3  was  flown  in  the  warm  sector 
of  a  cyclonic  storm.  The  soundings  obtained  from  the 
aircraft  data  indicated  that  the  atmosphere  was 
slightly  more  stable  than  moist  adiabatic  and  was 
saturated  in  all  cases.  The  height  and  temperature 
ranges  for  the  three  spirals  are  given  in  Table  1.  No 
significant  supercooled  water  was  detected  in  any  of 
the  storms. 


b.  Loop-averaged  spectra 

Snow-size  spectra  averaged  over  a  complete  loop 
are  shown  in  Fig.  4  for  the  first  spiral.  The  spectra 
are  labeled  consecutively  from  the  top  to  the  bottom 
loops  of  the  spiral.  The  spectra  are  approximately 
exponential  in  form  such  that 

j\(D)  =  JV0e_w>,  (8) 

where  N{D)AD  is  the  concentration  in  the  diameter 
interval  [D,  D  +  AD],  /V0  is  the  intercept  and  X  the 
distribution  slope.  Here  D  is  the  actual  particle  size 
measured  by  the  laser  imaging  probe. 

The  spectra  from  the  first  spiral  show  an  initial 
monotonic  increase  in  the  intercept  while  the  slope 
remains  relatively  constant  through  loop  12,  after 
which  there  is  a  rapid  decrease  in  both  the  slope  and 
intercept  during  loops  13  and  14.  The  remainder  of 
the  loops  show  essentially  no  change.  The  vertical 
separation  between  successive  loops  is  200  m. 

Because  of  the  quasi-erponential  behavior  it  is 
convenient  to  characterize  the  spectra  by  N0  and 
X  which  can  be  done  via  a  least-squares  fit.  Fig.  5 
shows  data  from  the  three  spirals  in  log,oN0-logioX 
(hereafter  N0-\  space).  The  loop  numbers  are  in¬ 
dicated  next  to  the  data  points  which  are  connected 
sequentially  by  a  line.  Environmental  temperatures 
are  indicated  for  places  where  significant  changes 
occur. 


Fio.  4.  Snow  size  spectra  for  spiral  I  (25  February,  1980)  averaged  over  the  various  loops 
Numbers  on  the  right  indicate  loop  numbers. 


178 


702  JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES  Volume  39 


•O  OO 

*  lew'1 1 

Flo.  5.  Nr  A  trajectories  for  the  three  spirals.  See  text  for  details. 


This  type  of  spectral  display  is  convenient  because 
a  straight  line  in  JV0-\  space  corresponds  to  a  con¬ 
stant  moment  of  an  exponential  distribution,  since 
the  yth  moment  of  a  spectrum  is 

Afj  =  D‘N#DdD  -  •  (9) 

Hence  in  logarithmic  N0-\  space,  the  moment  M, 
is  constant  along  any  straight  line  having  slope  ( j 
+  1). 

The  behavior  of  the  spectra  in  Fig.  4  can  be  dis¬ 
cussed  in  terms  of  the  trajectory  of  the  spectral  evo¬ 
lution  in  JV0-A  space  in  Fig.  5.  The  first  spiral  is 
characterized  by  a  gradual  increase  in  /V0  and  a  slight 
decrease  in  A  such  that  the  slope  of  the  /V0-A  tra¬ 
jectory  is  negative,  implying  that  all  spectral  mo¬ 
ments  are  increasing.  At  loop  12  the  spectral  evo¬ 
lution  changes  dramatically  and  both  N0  and  X  de¬ 
crease.  Eventually  N0  and  X  assume  approximately 
constant  values.  The  second  spiral  shows  a  similar 
pattern  of  three  stages  of  evolution.  Stage  1  is  char¬ 


acterized  by  an  increase  in  N0  accompanied  by  rel¬ 
atively  little  change  in  X.  Stage  2  is  characterized 
by  a  rapid  decrease  in  both  N0  and  X.  Stage  3  is 
marked  by  an  apparent  cessation  of  spectral  evolu¬ 
tion.  The  third  spiral  only  reveals  stages  2  and  3. 

The  JVo-A  trajectories  for  the  second  stage  of 
growth  are  roughly  parallel  for  all  three  spirals,  hav¬ 
ing  slopes  ranging  from  1.80  to  1.9S.  This  suggests 
that  during  this  phase  of  spectral  evolution,  the  sum 
of  the  diameters  of  snow  particles  is  a  conservative 
property  of  the  distribution.  These  three  cases  also 
suggest  that  the  distribution  slope,  X,  has  a  minimum 
value  of  ~I0  cm"1  which  characterizes  the  third 
stage  of  evolution. 

c.  Horizontal  gradients  in  the  microphysical  struc¬ 
ture 

In  order  to  examine  the  horizontal  gradients  in  the 
size  distribution,  the  mean  diameter  of  all  particles 
measured  by  the  200-Y  spectra  as  a  function  of  air¬ 
craft  magnetic  heading  is  plotted  in  Fig.  6  for  the 
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Loop* 


Fig.  6.  Average  particle  diameter  as  a  function  of  magnetic 
heading  for  spirals  l  and  2.  Diameters  in  mm. 


1st  and  2nd  spirals.  Each  loop  is  represented  by  a 
separate  graph  and  the  loops  are  stacked  vertically 
in  accordance  with  their  height.  The  approximate 
diameter  of  each  loop  is  6  km.  The  loop  in  which  the 
spectra  transformed  from  first  to  second  stage  growth 
(the  peak  in  the  N0-X  trajectory)  is  indicated  by  a 
star  in  each  case.  Before  the  transition,  the  mean 
diameter  is  essentially  uniform  within  each  loop  and 
gradually  increases  with  depth.  (Not  all  upper-level 
loops  are  shown).  However,  after  the  transition,  hor¬ 
izontal  inhomogeneities  develop  very  rapidly.  Note 
that  the  features  are  correlated  from  one  loop  to  the 
next. 

In  order  to  examine  the  point-to-point  behavior 
of  the  spectra,  we  arbitrarily  divide  each  loop  into 
four  quadrants  bounded  by  the  cardinal  directions. 
Ideally  the  spectral  evolution  can  be  studied  by  ex¬ 
amining  the  height  evolution  of  spectra  averaged  over 
a  particular  quadrant  rather  than  over  an  entire  loop. 
Fig.  7  shows  an  example  of  the  /V„-A  evolution  for 
the  second  spiral  for  the  four  quadrants.  The  fV0  scale 


is  different  for  each  quadrant  in  order  to  separate 
the  four  quadrants.  Note  that  the  four  N0-\  trajec¬ 
tories  all  show  the  same  general  features.  However, 
the  west-north  quadrant  starts  rapid  stage  2  growth 
earlier  than  the  other  quadrants  (e.g.,  examine  loops 
9  and  10  in  Fig.  6). 

Fig.  7  illustrates  that  even  if  we  examine  portions 
of  a  loop  the  spectral  evolution  is  coherent.  This  is 
consistent  with  the  previous  discussions.  Also,  the 
rapid  development  of  horizontal  gradients  of  mean 
diameter  are  apparently  related  to  the  fact  that  spec¬ 
tra  in  different  regions  undergo  the  transition  from 
stage  1  to  stage  2  at  different  heights. 

5.  Physical  interpretation  of  ASD  profiles 

The  regular  behavior  of  the  size  spectra  suggests 
that  the  ASD  technique  is  sampling  a  well-behaved 
physical  process  in  a  coherent  manner.  In  this  section 
we  will  discuss  the  physical  implications  of  the  ob¬ 
servations. 
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Fig.  7.  yr\  trajectories  for  the  four  cardinal  quadrants  for 
spiral  2.  The  Na  scale  refers  to  the  north-east  quadrant.  The  other 
quadrants  are  shifted  downward  by  one-half  decade  each  to  dif¬ 
ferentiate  the  data. 

The  effect  of  deposition  growth  on  the  size  distri¬ 
bution  is  characterized  by  how  the  rate  of  change 
of  the  particle  diameter  depends  on  the  particle  di¬ 
ameter,  i.e., 

fee*.  (.0) 

For  6  >  0,  larger  particles  grow  more  quickly  in 
diameter  space  and  one  would  expect  the  distribution 
slope  to  decrease.  For  &  <  0,  deposition  growth  will 
cause  the  distribution  to  steepen.  For  snow,  the  par¬ 
ticle  mass  varies  approximately  with  the  square  of 
the  particle  diameter  and  the  rate  of  change  of  mass 
is  directly  proportional  to  the  particle  diameter  so 
that  &  *  0.  To  a  first  approximation,  all  sizes  will 
grow  at  the  same  rate  in  diameter  space  and  the 
distribution  slope  will  remain  constant.  The  intercept 
will  increase  since  the  smaller,  more  numerous  par¬ 
ticles  will  grow  into  larger  sizes.  In  the  absence  of 


a  source  of  ^mall  particles,  a  lower  limit  to  the  size 
distribution  will  develop. 

The  first  stage  of  particle  growth  (observed  for 
spirals  1  and  2)  is  consistent  with  deposition  growth. 
This  is  not  to  say  that  aggregation  is  not  occurring, 
but  that  deposition  is  dominating.  Note  that  during 
the  first  stage,  the  increase  in  N0  and  the  relatively 
slow  decrease  in  X  suggests  that  the  concentration 
increases.  This  is  probably  due  to  the  growth  of  small 
particles  into  detectable  sizes.  The  source  of  these 
small  particles  is  either  nucleation  and/or  secondary 
production. 

The  second  stage  of  growth  (in  which  both  N0  and 
X  decrease)  was  observed  for  all  three  spirals.  This 
is  characteristic  of  aggregation  which  depletes  small 
particles  and  creates  large  ones.  The  ASD  data  for 
spirals  1  and  2  show  a  very  sudden  transition  from 
the  dominance  of  deposition  to  aggregation.  During 
the  deposition  growth  phase  the  particles  are  small 
and  collisions  are  rare.  Eventually  deposition  pro¬ 
duces  a  sufficient  number  of  large  particles  for  ag¬ 
gregation  to  commence.  Once  started,  the  large  par¬ 
ticles  produced  by  aggregation  accelerate  the 
aggregation  process  and  rapidly  deplete  the  smaller 
particles.  This  accounts  for  the  sudden  transition  and 
the  very  rapid  evolution  of  the  size  distribution  after 
the  transition.  For  spirals  1  and  2  the  transition  oc¬ 
curs  at  -15°C  which  is  in  the  dendritic  growth  tem¬ 
perature  range.  This  implies  that  the  tendency  for 
dendrites  to  form  aggregates  (e.g.,  Jiusto  and  Wcick- 
mann,  1973)  may  play  a  role  in  the  transition  to 
aggregation  growth.  However,  the  model  calcula¬ 
tions  described  by  Passarelli  (1978a,b,c)  do  not  re¬ 
quire  a  change  in  the  particle  geometry  or  collision 
efficiency  to  simulate  the  first  two  stages  of  growth. 

The  most  puzzling  behavior  is  the  apparent  sudden 
end  to  the  decrease  of  N0  and  X.  The  spectra  cease 
evolution  when  the  slope  reaches  ~10  cm'1.  This 
suggests  that  the  depletion  of  small  particles  is  bal¬ 
anced  by  a  production  mechanism,  and  that  large 
particles  produced  by  aggregation  are  somehow  de¬ 
pleted.  This  is  consistent  with  particle  breakup.  An¬ 
other  possibility  is  that  the  cessation  of  evolution  is 
merely  an  instrumental  artifact,  but  we  are  at  a  loss 
to  explain  how  this  could  account  for  the  observation 
that  all  three  spirals  evolve  to  the  same  limiting  slope. 

The  fact  that  all  the  spirals  evolve  to  the  same 
slope  suggests  a  collisional  breakup  mechanism.  This 
hypothesis  is  drawn  from  previous  work  on  drop  co¬ 
alescence  and  breakup  which  demonstrates  that  co¬ 
alescence  and  collisional  breakup  lead  to  equilibrium 
size  distributions  that  have  the  same  slope,  regardless 
of  the  precipitation  rate  (Gillespie  and  List,  1976; 
Srivastava,  1978). 

Recent  studies  of  ice  particle  breakup  (Hobbs  and 
Farber,  1972;  Vardiman,  1978)  emphasize  that  the 
coexistence  of  dense,  rapidly  falling  graupel  along 
with  fragile  aggregates  and  crystals  may  be  neces- 
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sary  for  collisional  breakup.  However,  our  observa¬ 
tions  suggest  that  the  presence  of  graupel  is  not  re¬ 
quired  for  breakup. 

Another  feature  of  the  data  which  could  be  ex¬ 
plained  by  breakup  is  the  tendency  for  the  small  par¬ 
ticle  concentration  to  increase  downwards.  This  was 
detected  by  a  1 -dimensional  imaging  probe  for  the 
small  sizes  (20-300  nm).  The  downward  increase  in 
the  particle  concentration  occurred  for  spirals  1  and 
2  and  is  surprising  since  aggregation  and  deposition 
both  deplete  small  particles.  In  the  absence  of  su¬ 
percooled  liquid  water  and  at  higher  temperatures 
at  the  low  levels  it  is  not  likely  that  nucleation  could 
have  been  producing  the  small  particles  in  concen¬ 
trations  of  order  10  L"1.  Herzegh  and  Hobbs  (1980) 
reached  the  same  conclusion  based  on  similar  ob¬ 
servations  in  similar  storms,  i.e.,  breakup  must  be 
the  source  of  these  numerous  small  ice  particles. 

A  detailed  quantitative  comparison  between  the 
observed  spectral  evolution  and  model  results  will  be 
presented  in  a  later  paper.  It  is  worthwhile  to  mention 
that  the  first  two  stages  of  growth  are  depicted  in 
Passarelli’s  (1978a,b,c)  analytical  model  of  deposi¬ 
tion  and  aggregation.  Incorporating  breakup  into  this 
model  is  not  difficult,  but  neither  the  probability 
of  breakup  nor  the  size  distribution  of  fragments  are 
known  for  snow.  Our  modeling  results  thus  far  show 
that  one  can  select  a  breakup  probability  function 
that  increases  with  particle  size,  and  be  assured  of 
reaching  an  aggregation-breakup  equilibrium.  How¬ 
ever,  it  seems  unlikely  that  one  can  use  the  observed 
spectral  evolution  to  determine  the  probability  of 
breakup  and  the  size  distribution  of  fragments. 

6.  Related  observational  studies 

Some  recent  observational  studies  of  snow-size 
spectra  have  shown  a  spectral  behavior  similar  to 
the  ASD  results,  although  the  scatter  is  typically 
much  greater.  Passarelli  ( 1978a, c)  employed  similar 
instrumentation  but  a  different  vertical  sampling 
scheme.  The  aircraft  was  placed  in  either  a  constant 
ascent  or  descent  while  flying  a  constant  heading. 
Particle  size  distributions  were  averaged  over  15  s 
intervals.  The  spectra  from  flights  on  6  and  10  March 
1975  are  shown  in  /V0-A  space  in  Fig.  8.  While  spec¬ 
tra  obtained  via  the  ASD  technique  show  a  very  sys¬ 
tematic  behavior  with  height,  the  spectra  obtained 
on  these  two  days  do  not.  The  spectra  apparently  lie 
on  a  line  corresponding  to  the  second  stage,  but  the 
position  is  random.  The  spectra  for  each  day  arc  very 
well  differentiated  in  N0-X  space  which  probably  re¬ 
flects  the  very  different  environmental  conditions  on 
the  two  days. 

Passarelli  ( 1978c)  also  averaged  particle  size  spec¬ 
tra  over  10  km  horizontal  passes  at  various  altitudes. 
The  data  for  a  flight  on  26  November  1975  are  also 
shown  in  Fig.  8.  The  passes  were  spaced  at  600  m 
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•  6  Mar  1975  tPo**ar*Ui  (978a  I 

•  10  Mar  1975  (  "  "  > 

*— i  26  Nov  1975  (  ••  1978(1 ) 

A  22  Jon  1976  (Houz*  et  at,  1980) 

□  8  Dec  1976  (  «  «l  al  1980 ) 

FlO.  S.  Related  observational  studies.  See  text  Tor  details. 

and  timed  at  10  min  apart — approximately  the  time 
required  for  snow  to  fall  from  one  level  to  the  next. 
The  results  show  a  first  and  second  stage.  Vertical 
incidence  radar  measurements  on  this  day  indicated 
fairly  steady  precipitation  while  the  aircraft  was 
sampling. 

Houze  el  al.  ( 1 979, 1 980)  employed  a  similar  tech¬ 
nique  by  flying  level  passes  although  no  attempt  was 
made  to  follow  snow  from  level-to-level.  The  spectra 
for  flights  on  22  January  1976  and  8  December  1976 
are  also  shown  in  Fig.  8.  The  data  are  more  scattered, 
perhaps  because  transient  environmental  conditions 
are  manifested  over  the  long  sampling  paths  that 
were  used.  However,  the  observed  spectra  are  still 
within  the  values  from  the  other  cases.  Although  the 
height  for  each  spectrum  is  not  indicated,  House  et 
al.  ( i979)  show  that  N0  and  A  decrease  with  increas¬ 
ing  temperature.  Based  on  the  general  trend  that 
temperature  increases  with  decreasing  altitude,  it  can 
be  inferred  that  /V0  and  A  decrease  with  decreasing 
altitude,  which  is  in  agreement  with  the  second  stage 
concept. 
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It  is  revealing  that  all  observed  values  indicate  that 
the  slope  X  is  always  >10  cm"1.  This  is  in  accord 
with  the  ASD  results  which  suggest  coliisional 
breakup.  No  matter  what  the  height  and  environ¬ 
mental  conditions  are,  the  observed  N0  values  are 
within  three  orders  of  magnitude  and  the  X  values 
are  within  one  order  of  magnitude. 

7.  Summary  and  concluding  remarks 

The  coherent  behavior  of  the  observed  size  spectra 
suggests  that  the  ASD  technique  can  successfully 
trace  a  1 -dimensional  region  in  which  the  effects  of 
spatial  and  temporal  gradients  are  minimized.  This 
technique  opens  the  door  for  more  detailed  quanti¬ 
tative  comparisons  between  microphysical  theory 
and  observation.  The  physical  interpretation  pre¬ 
sented  here  shows  that  snow  growth  evolves  through 
at  least  three  distinct  stages  which  we  have  chosen 
to  name  the  deposition,  aggregation  and  breakup 
stages.  Initially  deposition  dominates.  When  a  suf¬ 
ficient  number  of  large  particles  has  been  generated 
aggregation  becomes  dominant  and  produces  rapid 
changes  in  the  size  distribution.  Breakup  eventually 
limits  the  production  of  large  aggregates. 

The  observed  spectral  evolution  has  some  inter¬ 
esting  radar  implications.  The  transition  from  de¬ 
position  to  aggregation  might  be  detectable  by  ver¬ 
tical  incidence  radar.  The  transition  should  be  marked 
by  an  increase  in  the  height  derivatives  of  the  radar 
reflectivity  and  mean  fallspeed.  Deposition  will  gen¬ 
erate  a  fairly  gradual  increase  in  these  quantities, 
while  during  the  aggregational  stage  there  should  be 
a  rapid  increase.  Also  the  deposition-aggregation 
transition  should  magnify  horizontal  gradients  of 
radar  reflectivity  and  these  may  be  detectable. 

Another  interesting  aspect  of  the  observations  is 
the  rapid  aggregational  growth  which  should  produce 
changes  in  the  mean  particle  fallspeed.  If  it  were 
possible  to  perturb  the  deposition-aggregation  tran¬ 
sition  by  injecting  artificially  nucleated  ice  crystals, 
then  it  might  be  possible  to  exploit  this  for  the  re¬ 
distribution  of  snowfall.  The  ASD  technique  could 
be  used  to  check  this  hypothesis. 

There  is  a  definite  need  for  a  larger  data  base  and 
laboratory  studies  of  ice  particle  breakup  processes. 
Since  there  are  a  number  of  aircraft  equipped  to 


undertake  these  measurements,  it  is  hoped  that  a 
climatology  of  spectral  behavior  can  be  obtained 
through  the  efforts  of  independent  investigators. 
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APPENDIX  2 


LABORATORY  INVESTIGATION  QF  COLLISIONAL  BREAKUP 


In  order  to  investigate  how  the  collisional  fragment 
sizes  distribute*  a  laboratory  experiment  is  performed. 
Fig.  A.  2.1  shows  the  experimental  set-up.  The  main  apparatus 
is  a  household  freezer.  A  beaker  of  water  with  a  submerged 
heating  coil  serves  as  the  moisture  supply.  Two  pieces  of 
wire  are  mounted  on  clamps  and  are  placed  above  the  moisture 
supply.  One  of  the  clamps  is  of  a  movable  type.  Because  the 
wires  are  cold,  frost  crystals  will  grow  on  the  wire.  When 
the  crystals  are  large  enough  (at  least  3  mm)  the  moisture 
supply  is  removed.  A  Kodak  7302  blue-sensitive  fine  grain 
film  is  then  placed  underneath  the  wires.  The  movable  clamp 
is  then  gently  moved  so  that  the  crystals  on  it  will  collide 
with  the  crystals  on  the  fixed  clamp*  and  the  fragments  will 
fall  onto  the  film.  The  impact  speed  between  crystals  is  no 
mare  than  2  cm/s.  The  clamps  are  then  removed  from  the 
freezer  and  a  blue  strobe  light  is  flashed  directly  above  the 


film 

(Edgerton. 

1981). 

The 

crystal 

fragments 
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film 

cast 

shadow 

on  the 

film 

which 

can 

be 

measured 

after 

development  and 

enlargement. 

Since 

fine 

grain 

film 

is 

used* 

185 


enlargements  up  to  ton  times  con  bo  motfo.  Examples  of  these 
are  shown  in  Figs.A.  2.  2  ond  A.  2.  3. 

This  experimental  procedure  he«  boon  repeated  many 
timet.  The  results  from  typical  breakup  events  ore  shown  in 
Tables  A.  2.  1.  It  shows  that  the  number  of  fragments  is  high 
and  the  fragment  distribution  is  somewhat  negative 
exponential.  Therefore*  the  breakup  formulation* 
S(  x  l  x,  *  xJL)dx  ■  (x,  +Xj)  A’e'^  dx*  is  reasonable. 

This  experiment  has  the  drawback  that  crystals  are 
forced  to  break  up  and  no  aggregation  is  allowed.  In  this 
respect*  this  makes  the  experiment  not  too  realistic. 
However*  this  simple  set-up  does  provide  valuable  information 
on  the  relative  breakup  size  distributions  and  so  it  has 
served  its  purpose.  The  many  small  fragments  resulting  from 
collisional  breakup  are  especially  significant  in  cloud 
physics.  These  fragments  can  serve  as  centers  for  crystal 
growth  and  can  also  act  as  agents  in  the  generation  and 
distribution  of  electric  charges  (e. g.  Latham*  1963i  Latham 


and  Stow.  1967). 


Maximum  one-dimensional 

Number  of 

fragments 

fragment  size  (micron) 

Case  1 

Case  2 

0  -  200 

28 

12 

200  -  400 

14 

7 

400  -  60  0 

9 

7 

600  -  800 

3 

4 

800  -  1000 

3 

4 

1000  -  1200 

2 

2 

1200  -  1400 

0 

1 

1400  -  1600 

1 

2 

1600  -  1800 

2 

3 

1800  -  2000 

1 

2 

2000  -  2200 

0 

1 

2200  -  2400 

1 

0 

2400  -  2600 

1 

2600  -  2800 

0 

0 

2800  -  3000 

1 

l 

3000  -  3200 

(;  1  I 

i 

3200  -  3400 

2 

i 

3400  -  3600 

1 

0 

3600  -  3800 

1 

0 

3800  -  4000 

4000  < 

1 

3 

1 

1 

Table  A.2.1:  Samples  of  fragment  size  distribution 


resulting  from  laboratory  breakup 
experiments 
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